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I. INTRODUCTION 


The availability of @lebal Satel|]ite one lourag isms 
daily basis has exposed meteorologists to Weathened is eign. 
ances which were previously not known to exist. In 


particular, systems Called “polar lows (Hamleveu goo, 
nave been found to exist primarily over the oceans inthe 
winter seasons (Reed, 1979; Rasmussen, 1979). These storms, 
which have horizontal scales of Only 500— US 00 ieee cameo 
quite violent and occasionally exnibit thtnder Scene ee 
hail and hurricane-force wind gusts (Stevenson, 1968; 
Suttie;, 1970; yal oe 

[It has been suggested that two Gistincelwed) teoneme 
types of polar lows exist (Rasmussen, 1981; Locatelli 
at al., 1982; Sardie and Warner, 1983). The first of 
these is the Atlantic polar law yy when t yo ically sec sees 
Spiral structure and occasionally ane eve > hic Oo ces 
found in tropical storms. [iis edi stipe Cas spose 
norizontal scale of 500-1000 km and maximum amplitude at 
the surface. In addition, 1& tends to) farm wit i gece orn 
air mass north of the polar jet in an area with weax or 
negligible baroclinity. The second type 1s the Paciine 
polar low. This disturbance exhibits wavelengths of 1000- 
1500 km and is often a deep system, occasionally exhibiting 


height perturbations at 300 mb larger than those at the 


ik 


eee need, 1979) Sine Pacific polar low tends to form 
foeeenme Or tne polar jet in a strong baroclinic zone within 
a polar air mass. The overall appearance of this system 
mesemeres that Of typical Synoptic-scale extratropical 
Syseems, but of smaller scale. 

The. purpose of this study is to investigate the 
foaoenoe 10m the Growth of the Pacific polar low. It has 
been proposed by Reed (1979) and Mullen (1979) that either 
iemeoeropic instability, Conditional Instability of the 
Scomae<ind (CISKi@ur baroclinicw@instability may be respon- 
mee either singly or im some combination. Although the 
conditions for barotropic instability are often met in the 
aoeireeatring polar low formation (Mullen, 1979), Gall 
(1976a) showed that barotropic effects are about an order of 
magnitude smaller than baroclinic effects for wavenumber 15. 

Many analytic studies support the argument that 
Sees linive instability is responsible for the growth of 
the Pacific polar low. The earliest studies of baroclinic 
iwedb lity “O@rarney, gay: Eady, 1949) indicate that 
4000 km waves are most unstable and observations verify tnat 
Giiseseale GdOmifates in extratropical regions. However, 
iementeG@ent analytic studies snow that the growth rate 
spectrum is very sensitive to the Richardson number (er) e 
immetotseene ratio of tne Static stability to ghe square of 


the vertical wind shear. Reduced dete tee lowest levels 
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tends to shift the maximum growth rate to shorter wave- 
lengths (Gambdo, 19/70a; Tokioka, 19702 Mansi 1] yam? eee 
1976a; Duncan, 1977; Staley and Gall, 1977). However, the 
role of baroclinic instability 1s Still Wn “gomibt “s ieee: 
baroclinic waves tend to be quite shallow, as opposed to tne 
‘observed deep structure of the Pacific pole! oW- et ecg 
tion, Sardie and Warner (1983) used a quasi-geostrophic 
model to show that the maximum growth rate for conditions 
typical of those during the fommetioneoh Sameme cc ii Gaels 
low 18 fOuUnd at @ Wavelenugtn Gia 00 in 

The presence of convection suggests that it may be a 
Significant factor in the development of the Pacific polar 
low. Rasmussen (19/79) s has Shiown se iidtest [SK sis camden 
producing "“mini-hurricanes" similar to the Atlantic polar 
ike We “cecal CISK requ res.daspreseoas timiG, ¢ 1 relation 
whicn must be explained. Also, CISK requires a small 
vertical! wind shear, a Condition ipso) cen tie ee npiCe ea 
development of the Pacific polar low (Mullen, 1979). 

Therefore, it appears that some combination of 
baroclinic..instabilityeand didibathiGmarocesicics ise ees 
convection, are responsible for thesdeye lope emote: 
Pacific polan lows It Nas@bielen@siiloiy Tee ielesRe mco lias iiariimrd 
CISK-type latent heating shifts the most unstable wavelengtn 
from the synoptic-scale to 1000-2000 km waves (Gambo, 1976; 


Rao and Ferreari, 1979; Sardie and Warner, 1983). Also, 
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Gall (1976c) demonstrated that convection can serve to 
Mileeed>encne Vertiedl extent of medium-scale (2000 km) 
waves. 

Sigenace flUxes Of heat and moisture also must be 
considered Since large air-sea temperature differences are 
SeeemedssOciated with Pacific polar lows (Mullen, 1979). 
Large upward fluxes of neat can reduce the low-level static 
Stability which, as explained above, can result in baro- 
clinically unstable 1000 km waves. The role of surface 
muteces in Strengthening disturbances has been discussed by 
Petterssen (1956), Danard and Rao (1972) and Forbes and 
eie@iete, 1984). 

Based on the observational and analytic studies 
described above, the major hypothesis of this study is tnat 
Pacific polar lows are formed primarily through baroclinic 
instability. They are a preferred mode, especially in 
winter, because of the unique vertical temperature profile 
forced by a polar air mass moving over relatively warm 
water. The upward sensible heat flux in the polar air 
reduces the low-level static stability and therefore, Re. 
This iS precisely what is needed to set up baroclinically 
unstable waves on the scale of the Pacific polar low 
according to the theoretical studies discussed above. It 
Pomccommmo thesized thas surrace fluxes and latent heat 
release through convection modify the baroclinic wave Dy 


Bisener IN tensitying it and increasing 16s vertical extent. 


21 


To test this hypotnesis, the numerical wWo@e)) Usman 
Fleet Numerical Weather Center's Navy Operational Regional 
Atmospheric Prediction System (NORAPS) is used (Appendix A). 
This model includes convective and non-CcOnNvective preci ime. 
tion processes and state-of-the-art planetary boundary layer 
(PBL) physics. In addition, the model wWdy spemgun wie 
Various combinations Of "Grid soacinGg sande. cat) or en nee 
makes the model quite useful for studying waves of various 
Sizese 

Interpretation of model wesults is decomp lysed ear, ana 
budgets of kinetic and available potential energy as given 
in Robertson and Smith (1983). These equations (Appendix 8B) 
are split into zonal and eddy components and the exact form 
of the available potential energy 1s Wisec. Sites tie mnomes 
cast model uses a closed domain, the eddy components can be 
directly related” to a Single WwaVe Im sume forecedst scamaim- 

This study consists Of twowmage@r SeCCtlOmsS.  nmemen | 7 sce 
Chapter [I1, describes linear growtherates Comoe eas ice. 
linearized eh wes the NORAPS model. Two sets of initial 
conditions are used. The first is a mean winter jet and tne 
second is a jet typical of conditions during the development 
of Pacific polar lows: Small SeneUreat Onseare Inerganc ed 
into the initial fields and an ddiapatre= vegs (oso pene 
model is integrated while restoring the initial mean zonal 
fields every time step. The effects of linearized friction 


and wave-zonal flow interaction darew discuss ecasuinitomrs cat gre) 
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aisOesnows NOW Ddrotropie and Daroclinic instabilities 
contribute to the generation and destruction of the eddy 
Kinetic energy of tne most unstable waves. 

Ime second major section, Chapter IV, describes 
nonlinear model simulations of the polar low. Detailed 
analyses of wave structure and energy budgets are presented 
Homes abatic and didabatic simulations. Tne sensitivity to 


tne air-sea temperature difference is discussed. 


ll HUMERT CARP ROCEDURiD 


The NORAPS forecast model has several features which 
make it applicable for extensive numerical simulations. 
First, the domain size and the horizontal and vertical 
resolution in the model may be varied. Second, dynamics and 
ohysics appropriate to subsynoptic fedtures dre ime miaed sa. 
described in Appendix A. Exact values for various param- 
eters are given in the description of each model experiment. 
A map factor of 1 1s uSed OVer Che engl reed oman cece 
forecasts. Also the Coriolis force, Tf, 1s “eCaleutaveemisan, 


the p-plane approximation, 


aig. Oe) (22s 


A. INITIAL DATA 
eee odbc elie ice ves 

Two baroclinicdally unstable jet profiles werevenosce 
TO eis Gi uicly. The first jet describes the mean winter 
conditions over the Pacific. The conditions were derived by 
Lau (1978), who composited National Meteorological Center 
(NMC) analyses from 850 £0 MOQ Wo tone wig oe nee 
sector from 150E to 150W and from 20N te JON sia ene oie 
included 1000 mb winds estimated by Crutcher and Meserve 
(1970). It was expected that this jet (referred to here- 
after as the mean winter jet) would be appropriate for 


Studying typical oceanic synoptic-scdaie lows. 
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ite second jeteiaemerred to as the genesis jet, is 
typical of conditions during tne genesis period of north- 
eaoberm Paciiic polar lows (Mullen, 1979). Mullen (1979) 
Tenimedetnis period to be “at the latest hour, either 0000 
GMT or 1200 GMT, before the cloud mass in question showed 
any discernable comma or hook shape in the infrared (IR) 
Pagery., Mullen composited NMC analyses for 22 north- 
masotern Pacific polar lows. In all cases, the width of 
mieeclOUd mass at the genesis period was at least 500 km. 
It was expected that this jet would be appropriate for 
seuaying the polar low. 

The mean winter and genesis jets of Lau (1978) and 
Mullen (1979) were transformed to model coordinates in the 
following manner. First, the channel width of the model was 
emoscem te pe 6400 km with the jet core in tne center of this 
eianm@e!. Second, at all Cero ove lomo EO 00 To aN 
the genesis jet and up to 100 mb in the mean winter jet, 
wind speeds were tabulated in 250 km increments over an 
fieerior.cs00 km meridiomal slice. It was assumed that the 
Boaiewind £Or DOLN jets decreased linearly to zero at both 
the north and south walls. The shape of the genesis jet 
woe somo Was assumed Lo resemble that of the mean winter 
jet. Third, these wind speeds were linearly interpolated 
Dee cenwemandaqcOnry pressure surface Co the model horizontal 
Tardepormes. Fourth, Ome pass WIth “aneelepotnt smoothner- 


desmoother (Shapiro, 1970) was applied to these fields to 


Ze 


eliminate any small perturbations which may have been 
introduced in tne tabUuUladeion or Interpo las onme | sc nomena 
Cross-sections of the mean winter and genesis jets are shown 
| rc 2 ele 
Zen © del tezrarcaenn 

Tne equations used in the initialization are greatly 
Simplified since initially, all x-derivatives are zero and 
v=Q over the entire domain. The initial north-south surface 
pressure gradient is computed geostrophically from the v- 


equation Of -mMoOUtTONe im a Goo nd Ia uc oe 


d91nt = fu (2.23) 





subject to the boundary condition that the Gradaaen) 07 era 
zero across the boundaries. A constant surface temperature 
of 282 K is used and the 1000 mbeu-components sine; 1o3e2od 
are specified. The APE mean surface pressure is specified 
as 1013.25 mb. Initial geopotential sweep res sim ous eres 


are computed geostrophically using 


cl ie | 
Sy fu Necc|y) 


Subject to the boundary condition thicttect hom qciucy) emoumecomes 
zero across the boundaries. The mean height of each level 
is computed as follows.  Veimpemaletie sence ees the mean 
winter and genesis cases, as given by Mullen (1979), are 
tabulated and interpolated to grid points in the manner 


described in the previousssect \om for, tne. zonal ewig 
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Heights are computed hydrostatically and the areal mean is 
found by averaging over each level. 

The initial fields are then interpolated linearly 
with respect to p0.2051 to the model o coordinates. When 
the model is integrated with these fields, a north-south 
pressure oscillation, or “sloshing” of 3-4 mb is observed. 
To achieve a better balance mmet ness) tae 1 clei cemeneny 
height fields are computed using the v-equation of motion 


in Go Coordinates. 


(2.4) 


CK en - 31nd 
3 fu-RT 3y 


Subject to the boundary condition that the raggeneno) o 1s 
zero across the bDOUmdayries. 8 far (2e4) mMWet be Sei odeay an 
iterative process, since the temperature appears on the 
rignt side. Qinieitemriepcremecisce the temperatures obtained 
from Eq. (2.3) are used. On each succeeding pass, tempera- 
tures obtained from the previous iteration are used. The 
procedure 31s repeated for flVe Wterdgy ons. ec seo eng 
the fields directly on the modelc coordinates, the pressure 
sloshing was reduced to less than 0.1 mb. After the 
Initialization, the genesis fields exhibit superadiabatic 
lapse rates since no formal Constraints: sung e so -aem On eee 
genesis case, the model is integrated for 12 h neglecting 


friction and diabatic effects to achieve a balanced state 


without superadiabatic layers present. Cross-sections of 
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the initialized temperatures for the mean winter and genesis 
eases are SnOWwn in Frq. 2-2. 

Picea Indl wonoemostlec Variable to’ be initialized is 
tne moisture. For both jet profiles, the zonal mean winter 
relative humidity fields for the northern hemisphere from 
London (1957) are used. Values are tabulated and interpo- 
lated in the same manner as for tne zonal wind. The 1000 mb 
relative humidity given by Mullen (1979) for tne genesis 
meage Of the polar Vows is used. A cross-section of the 
woteral relative humidity is shawn in Fig. 2-3. 

The remaining variable to be specified is the 
Sea-surface temperature (SST). All experiments apply for 
meteerater Conditions. Tne specific SST distribution is 
described with each experiment. There is no east-west 
Gradient in the SST and the values are fixed throughout the 
forecast. The surface moisture is the saturation value 


S@eeseponding to the SST. 


B. ENERGY BUDGET 

An energy budget is useful for studying tne development 
and maintenance of the atmospheric disturbance. It is 
common practice to decompose the energy equations into 
Saiaerons fOr zonel @nd eddy components of the flow. This 
mmpanreiculdrily useful in this study, as @ perturbation 
iewscein a zoOnd|l floweeEnegdy transtems within the eddy, 


and interactions between tne eddy and the zonal flow, will 
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be studied. The effects of friction and atebati Cee eee. 
also will be described in terms of the energy budgets. 

Many earlier studies have used the energy equations 
developed by Lorenz (1955), which include an approximate 
form of the available potential energy. Gall and Johnson 
(1971) demonstrated that this form of the equations can 
produce erroneous generation of eddy available potential 
energy. For this reason, an exact formulation of the energy 
equations, as derived by Robertson and Smith (1983), is used 
in the study. oA description of tne budget equations is 
Given in Appendix 8B. 

The energy budget terms are calculated using output 
fields from the forecast model. All fields are saved each 
6 h and are time-avereged during the model integral iaiger om 
a4 one-hour period centered on the oWeput Cime. (i ciieigeiaer 
a l2 > forecast field represents the average of all ‘forecast 
times between 11.5 and 12.5 hours. All fields dare imeer- 
polated linearly with respect to In p from the model g 
coordinates Mae. Constante pressure surfaces. There are 22 
constant pressure surfaces In SOSmMbe increments Starving edas: 
25 mb and extending down toe] 07 So) m DO eco sc Spl lentes meer 
or above the surface are used in the budget calculations. 
That 18, any mean Guat Gy iat.) 07 Secret Cragmcuianig 


only those points for which ¢ie VsiteGacempimesstre 7 omemioe 


a 


The transformation of the vertical velocity from the 
Mode! so coordinates to pressure Surfaces necessitated the 
saving of several fields. This transformation is obtained 
Pee osottueting EqQ.erA=l) Into @ = dp/de, which yields, 

beige © 6 2 aise | (2.5) 

at 

wiewitirst term om tie right side of Eq. (2.5) is simply 
the model vertical motion at the interfaces of the model 
meayers. 9 Ihese fields are interpolated directly to the 
Meeeoure surfaces. The last two terms are calculated at the 
middle of each model o layer and are interpolated directly 
irom these levels to the pressure surfaces. As a check, w 
was also calculated kinematically using the time-averaged 
momocievyerields on pressure Surfaces. Differences in the 
budget results were less than a few percent. 

The computation of the reference pressure is done as 
meron cem First, In p 1s interpolated linearly with respect 
meme urone une model ao coordinates to isentropic surfaces 
Palen tervalsSwe All 1Sentropic surfaces below the earth's 
Surface are assigned a value of In mw. Second, the areal 
feameeressiire (t.és, the reference pressure, oe is 
Semeeilatkea On Gach isentropic surface. Third, In Pp. las 
Piecmmoldled tO presstire SUrfaces linearly witn respect to 
peeemerT OSS-seclions Gueunaetn| tialmetTfi¢iency factor for both 


eeestriuctures are shown in FigGege2eee 
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No horizontal interpolation is necessary for the mode! 
Meow, fields. AN] horizontal derivatives are computed in a 
manner consistent with the model finite differencing on the 
soaiemomemgnid. Centeredsdifferencing is wsed for vertical 
derivatives. The time change in'energy is computed using a 


Gemcered difference over a 12-h interval. 
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ITT. tNEAR GROW iSite 


The concept of using the perturbation technique to study 
wave growth in a zonal flow was introduced by Charney (1947) 
and Eady (1949). These studies indicated that the most 
Mecanle wavelength for mean atmospheric conditions was 
about 4000 km. In general, this agrees with observations 
whicn show that the extratropical regions are dominated by 
waves of this size. However, more receme otudies mos o1 mom 
(1972) and Gall (1976a) indicate that the most unstable 
wavelength depends critically on the mean flow structure. Gall 
(1976a) showed with the linearized equations that wavenumber 
15 was most unstable for mean winter conditions. However, 
wavenumber 15 failed to dominate over longer wavelengths due 
to wave-zonal flow interaction dnd siieface friet jommen accuse 
for finite amplitude. 

The purpose of tnis chapter is to study the growth rates 
for the mean winter and genesis jets described in Chapter 
[Il. First, linear growth fates Wil | pemeomou tt denier both 
jets and the energetics and structure of the most unstable 
waves will be examined. Second, the effect of wave-zonal 
flow interact on wi ioe determined. silnmegdse theme (ec tac 
linearized surface friction on the growth rates will be 


examined. 
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Pree PERIMENTAL DESIGN 

A reduced version of the numerical mode! described in 
Appendix A is used to study the linear growth rates of 
Pemeun~Ddtk tons in the Mean winter dnd genesis jets. Two 
Simplifications are made to the model in these experiments. 
meeocy tie physical processes are not included. Second, 
after each time step, the zonal averages of u, v, T and nm 
are replaced with the initial zonal averages of their 
@eseective Values, which eliminates interactions between the 
Growing waves and the zonal flow. These See enone 
transform the nonlinear model into a linear model with 
@espect tO interaction witn the mean flow. The nonlinear 
interactions between waves are still included. 

Micmodtecsian grid Nas a tori zontal resolution of 100 kin 
with 94 grid intervals in the x-direction and 64 intervals 
in the y-direction. Therefore, the east-west extent of tne 
iercers 9400 km or 1/3 of the global circumference at 45N. 
Since the genesis jet exhibited the largest growth rates for 
Mmeeeciortest waves, anotmer version of the model was run for 
the genesis case only. This version uses 75 km resolution 
Pee  ifitervals in the x-direction and 85 intervals in tne 
Vememrection. Ihe @ast-west extent of this grid is 7050 km 
en 1/4 of the global circumference at 45N. The model 


vertical levels are given in Table A-l. 
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[Initial v-component pertunodat ons (aie Gp odie cas 
model wavenumbers 1-10. Tnese correspond to earth wave- 
numbers 3-30 (4-40) in increments of 3(4) fon tHe PU0Ry ean 
grid. The perturbation, whieh is constant in the vent aca 
is given by, 

yt =" “Sin BES) 54g? (2a |as) 

x 

The y-structure is chosen so v'=0 on the north and south 
boundaries. [ne Wil tigewieamel 1c demon ace eRew nz eelsneiherncs =| & 
OF sol for the mean winter jet. However, the initial 
amplitude for the genesis jet was set to 0.01 m s-1l because 
it exhibited such high growth rates. This insured a longer 
period.over which to eomoute. tie” lime cigs Gucci fei aiCi mers 
noted by Gall (19/764), the imposition of the initial epewoume. 
Dation leads to generation of gravity waves. Due to the 
presence of these waves, it was estimated tnat the growth 
rate was computed £0 Wi Giulia ede Glide clo: i el O) sunita ene 
Be ayalesy me S16 W|ta\ 7. 

The technique used to evaluate the linear growth rate is 
Similar to that used by Gall (1976b), who inserted only one 
wave per forecast. Therefore, Gall required many forecasts 
to achieve the growth rate spectrum for the given zonal 
flow. In this experiment, 10 wavenumbers are included 
initially, so that interactions between waves may alter tne 


growth rates of some of the waves. However, these nonlinear 
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Interactions will remain small since the perturbations are 
small. As a check, the model was run with only one wave- 
Pimiber in the inittal fields for two different wavenumbers. 
The Growth rates computed from these runs were identical to 
those computed when all 10 wavenumbers were included, which 
verified that the nonlinear effects were small. 

The growth rates and phase speeds were computed as 
follows. A double Fourier decomposition of the v-component 
ieeme WoOwest inodel layer was calculated every 12 h to get 
Picmampyitude and pmase of each of the 10 waves. During the 
period of linear growtn, the growth rate for each wave was 


computed as 


; ny Seley 
uC Se ee SS ee) 


1 at ANE 


where tne subscripts s and e refer to the start and the end 
of tne linear growth period. Phase speeds were also calcu- 
lated for each 12-h period, and were averaged during the 


period of linear growth for each wave. 


B. WAVE GROWTH 
oomedt Growum Rates and Wave Structure 
The growth rates computed for the mean winter and 
genesis jets are snown in Fig. 3-l. There are two major 
differences between the curves. First, the genesis jet 


Siro lus much larger Growth rates, especially for 


a8 


BY 


stsauab (q) pue wo qULM ued ( 


YIGWNN SAUM 
BE B2 BI 


=] 
e ) 


uauiLuadxa Y Yas so} sya 


AO} SOPRU YYMOUH sPauL] <q[-¢E 


YSEWNN SAM 


a-g | OY BE Bz 


cB 


y 


Oowvwrnrn BO DOD WO 


8° 
g° 


B 


—=—~ est = es RQ WB 
(ABO Ysd) SLBY HLMOYD 


J 
c 


BI 


40 


Wavenumbers greater than 20. Second, the most unstable 
wavelength (3 ManKeciemmarhfTerent fOr the two jeus. For the 
mean winter jet, wavenumber 12 is most unstable, with a 
eeovtnerdcve Of 0.59 day-l. For the genesis jet, wavenumber 
Piets most UnsStablewwith a growth rate of 1.82 day-!, and a 
secondary maximum of 1.29 day-! at wavenumber 8. Wave- 
numbers larger than 40, which were not geaced here, may 
Saoiteeven largermegrowth rates. 

The growth rates for the genesis jet were relatively 
easy to determine. However, growth rates for the mean 
winter jet for wavenumbers greater than 24 were more diffi- 
cult to estimate. These waves exhibited little or no growtn 
for up to 7 days and then experienced rapid growtn. The 
longer waves (wavenumbers < 24) generally began exponential 
Growth between days 2 and 3 of tne forecasts. It was judged 
that the eventual Growth of the shortest waves was due to 
Poencapinleraeraones with tne longer waves, which by day 7 
fmiredened amplitudess> 5 m s-1. Model runs including only 
wavenumber 28, or including only wavenumber 40 verified tnat 
the linear growth of these waves are quite small for the 
mean winter jet. 

Miewexplananion fOr thie Snorter wavelength disturd- 
cieewi tne genesis jet is illustrated yy een e Nl cumGis.OM 
number (R.) an S tigdiie 1 Oncmr One Dom ie etice (Ff agen 3-2 ) R ins 
miemiceto Of the Staeve Stab lityeto ene square of the 


vertical wind snear. Staley and Gall (1977) demonstrated 
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ieemesmat'l, Yow-}evel static stability and/or an increase in 
memo —-level wind shea (i.e., a small Rs) tends tomshift 
tne maximum growth rate to shorter wavelengths, It can be 
crim Omer 19. 3-2 thatthe mean winter case has Rs < 10 
Over a large depth Ct eeiesteopospiere in the Vicinity of the 
jet and in the low levels poleward of the jet. The genesis 
get case exhibits the same pattern We ee eae edie ee amOnW = 
imereeues poleward of the jet core. In fact, many values are 
teecmenan 1. This difference in R. 1S consistent witn tne 
large difference in growth rates between the mean winter and 
genesis jet cases for wavenumbers greater than 21. 

The normalized amplitude of the waves during the 
period of linear growth is similar for the mean winter and 
the genesis jet cases (Fig. 3-3). Wavenumber 6 has a maxi- 
mum amplitude in the lowest layer and also near 300 mb. The 
Shorter waves (wavenumber > 15) are confined to the lower 
troposphere with the maximum amplitude in the lowest layer. 
The amplitude of these short waves decreases to 20% of the 
meV eacs values ot 570 mb for theemean winter jet and at 
Sxoewmo tor the genesis jet. 

These results are consistent with those of Staley 
and Gall (1977), who used a 4-level quasi-geostrophic model 
momo bug ye tne erfect Of wind Shear a@mdestatic stabi ! ity on 
baroclinic development. They found that decreasing the low- 


PPestatyc  stdbi lity Hi-c., deCweacimask |eshifts the 
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wavelength of maximum growth to shorter wavelengths. The 
Geenerencese et ween tae growth rates presented here for 
BiesmedmM winter and genesis jets is evidence of this. 
Seavey wand Gall (1977) also folind that wavelengths less 
than 3000 km (wavenuniber >10) are confined below 500 mb 
and exhibit ener? maximum amplitude at tne surface. 

iim Biase = Strucmune Of eWavenummers 6, 12, lee 24 and 
30 for the mean winter jet and wavenumbers 6, 12, 18, 24, 30 
and 40 for the genesis jet are shown in. Figs. 3-4 and 3-5. 
For both jets, wavenumber 6 has a pronounced tilt to the 
west with height between tne surface and 700 mb and a 
Becucea westward tilt up t@ 300 mb. For wavenumbersei2 and 
18 in the mean winter case and wavenumber 12 in the genesis 
case, a moderate westward.tilt is observed - the lowese 
400 mb with little or no tilt -above. For the remaining 
shorter waves, a very strong westward tilt is found in the 
lowest 200-400 mb. This tilt becomes more pronounced for 
Increasing wavenumber. Since a westward tilt is associated 
with baroclinic processes, these results suggest that baro- 
GSrimicmemstability is exhibited by all] waves, but is 
confined to the lowest few hundred millibars for wavenumbers 
Siaeater Gnan 12. 

The phase speeds for the perturbations in tne mean 
winter and genesis jets are shown in Fig. 3-6. Phase speeds 
in the genesis jet case are 30-40% higher tnan those in tne 


mean winter jet as would be expected from the larger wind 
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Speeds in the genesis case. For the mean winter jet, the 
maximuin phase speed occurs at wavenumber 9, and there is 4 
fomerdicecredase in Phase speed for shorter waves. For the 
genesis Jet, the maximum phase speed is also at wavenumber 9 
with a similar decrease in phase speed for the shorter 
waves. The lower phase speeds for the shorter waves are 
because these shallow systems have a lower steering level. 
@. Energy Conversions 

The most unstable mode for the mean winter jet and 
the two most unstable modes for the genesis Jet are now 
examined in more detail. For each of these waves, the model 
momcomnigured to allow just that wave in the domain, so the 
process(es) leading to its development can be isolated. In 
particular, the baroclinic conversion (CE) between eddy 
available potential energy (AE) and eddy kinetic energy (KE) 
and the conversion (CK) between zonal kinetic energy (KZ) 
and KE are examined. 

The baroclinic term, CE, is the correlation between 
werewcal)s velocity and temperature perturbations, CRD Ese Ke 
mienedses fOr rising warm air and sinking cold air at the 
meme revation, and vice versa. Ihe conversion CK is 
Seaiadeed inte thlewhorizontal (HCK) and vertical (VCK) 
components. Simce v=0 in this linear study, HCK pepe scutes 
Ciempareoiropic term (io iin.) which is the correlation of 
u and v component perturbations multiplied by the meridional 


iememeeor the mean zonal wind. North of tne jet, where 
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su/ay < 0, a positive (negative) correlation between u' and 
v' leads to an increase ldecredse). in) Kb eee dilia. 
troughs which tilt opposite the horizontal wind snear are 


are barotropically unstable. VCK represents the term 





(-u'w!au/3p) and therefore represents the same type of 
motion as HCK but in the Vertical. 01 anc eee 3s and 
increases with height, a positive (negative) correlation 
between u- and e will lead to an increase (decrease) in KE. 
This is referred to as Kelvin-Helmholtz (shearing) instabil- 
ity which can become significant for ees (WE Oil ainis | 4 
1968). The term VCK is not included in quasi-geostropnic 
models such as the one used by Sardie and Warner (1983). 

The most unstadle mode for the mean winter jet is 
wavenumber 12, referred to as W12, wnich corresponds to a 
wavelength of about 2400 km. For the genesis jet, the most 
unstable waves are wavenumbers 8 (G8) and 40 (G40), which 
correspond to wavelengths of about 3500 and 720 km, 
respectively. Three different model domains are used to 
separately simulate these waves (Table 3-1). As in the 

Table 3-1. Number of grid intervals and horizontal 
resolution for waves W12, G8 and G40. 


Wave x-intervals y-intervals H Cigale Mitacal 
Resolution (km) 


W12 24 64 100 
G8 35 64 100 
G40 18 160 40 
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previous experiments, the physical processes are neglected 
meee ne zonal averages of u, Vv, 1, and w are restored each 
time step. 

[emt imevevolUcion Of the Lerms CEMEICI and Vick for 
W112, G8 and G40 are given in Figs. 3-7, 3-8 and 3-9, respec- 
mimeije FOr Wi2, CE is the dominant term, which shows that 
the growth of this wave is primarily through conventional 
Paggemminic iNStability. HCK and VCK are 5-6 times smaller 
than CE, which indicates that barotropic and vertical shear 
ecco cre small. Also, HCK and VCK tend to cancel since 
they are nearly tne same size but of opposite sign. For G8, 
CE iS again tne dominant term. VCK is positive and about 
One-third of CE, which suggests that Kelvin-Helmholtz 
instability may play a secondary role is the development of 
this disturbance. The HCK is again negative, but is only 
avoueemia!t the magnitude of VCK. For G40, both CE and VCK 
are positive and nearly equal in magnitude, which shows that 
this wave draws as much energy from KZ as from AE. Term HCK 
is negative and an order of magnitude less than CE and VCK. 

rt 1S apepamwemt from Pugs. 3-/.553-38 and 3-9 that the 
memaun Gf time of linear growth of G40 is significantly less 
than tnat of the longer waves G8 and W12. While the linear 
Growth rates extended past 120 h for G8 and W12, they only 
continued to about 60 h for G40. The reasons for the rapid 
Seedhization of this wave will be discussed in the 


mOlNNOWI Ng Section. 
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Fig. 3-7. Time variation of (GE, (SH@ andy @ eer nae 
number 12 in mean winter jet. See text for 
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mide o-oo. Simibar to Fig. 3-7 except for wavenumber 8 in 
genesis jet. 
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Fig. 3-9. Samal cigeneoue: icc except for wavenumber 40 
in genesis jet. 
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ioe oWeSt Cp OSs ese CL MOMomOn Veit led Wmevelocity, 
temperature and trougn position for W12, G8 and G40 are 
Caen fF 19S. 3-10, 3-11 and 3-12, respectively. Very 
Meas strucmunres dre Observed for W12 and G8. For both of 
these waves, tne trough tilts to the west with height with 
the maximum slope from 800 to 600 mb. The warmest air is 
moundedmedad Of the trough with a slight tilt to the east 
with height below 700 mb and a slight tilt to the west witn 
neight above 700 mb. The maximum upward motion is found at 
modeme-. A Droad area of relatively weak subsidence occurs 
behind the W1i2 and G8 troughs. 

IWe -structtire of the G40 wave is significantly 
Grtperent from Wl12 and G8 im that the warmest air slopes to 
tne west witn height, as does the trough. Therefore, only a 
emaerow layer (about 250 mb thick) of warm air is found 
above the position of the surface low with relatively cold 
air above 750 mb. The upward motion above the position of 
the surface low leads to the generation of KE below 750 mb 
mie ete > Orr and =a desimemcrmon of KE above /50°%!1b. tInis 15 
consistent with tne shallow nature of tnis wave. 

The relatively large values of VCK for G8 and espe- 
cially for G40 are due to the strong vertical wind snear in 
the slew aA Noes nce Go Cmand sl a0<Uemtmis implies a 


positive correlation exists between u' andw'. Fields of u' 
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ieee Poe oimilar te rig. “3=lO-except for wavenumber 8 
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ieee TOr G40 (Fig. 3-13) exhibit a pattern of wpward 
motion in areas of east winds and downward motion in areas 
of west winds. 

Orlanski (1968) showed that Kelvin-Helmholtz insta- 
cimney Cam dominate over conventional baroclinic instability 
Hor small R-. Stone (1966) showed that helical circulations 
Send ie| to the zonal flow dominate for 0.95 < Re < 0.29 
(symmetric instability), whereas the circulations are 
egeemebeuiar to the zonal flow for Ry < 0.25 (Kelvin- 
fone Zeanseadility). A north-south cross-section of v 
mierreennmougn G40 1S given in Fig. 3-14. Strong vertical 
snear in the v-component coupled with tne horizontal 
gradient of wleads to a circulation in the vertical plane 
about a horizontal axis near 800 mb. Since tnis is a north- 
mernememoss-Section, this circulation is parallel to the 
zonal flow and Suggests that symmetric instability is 
present. However, the u' and w' fields in G40 (Fig. 3-13) 
sgest thdt this disturOance 1S Strongly tilted by tne 
norizontal wind shear and therefore this helical circulation 
iaveaiso De driven by the Kelvin-Helmholtz instability. R 
values computed in the lowest level range from Q.03-2.0, 
Moe indicates tWat either type of motion 1S possible. At 
mivemdue, tne spresencewog ER 1S reldeal coer Cusleaukeon aes 
consistent with the results of Stone (1966) and Orlanski 
(1968). Therefore, based on the energy conversions and wave 


Beavceumens 1c 415 concluded that G40 is a “mixed” baroclinic 
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wave. Tnat is, energy is gained from both KZ and AE. 
Although VCK 1S Dario ReGiie "Ob atg eo 1c sees (eee? 
considered a baroclinic DroceSsmeqnce 1t imo lvesmer ie 


VErETed! wind sSilean:. 


C. EFFECT OF WAVE-ZONAM FEOW INTERACT DON 

Although Gall (1976a) showed that the most unstable 
Growth rate occurred at wavenumbers 12-15, the dominant 
scale of motion was eventually found to be wavenumbers 5-/. 
Subsequestly, Gall (1976b) included wav emzonial@in love Iie i= 
action and surface friction e@ffects.  Timawas fomopdetnde sen 
growth rate in the lower troposphere exhibited the most 
rapid reduction due €0 Gniereasimg Siatie Stability ana 
reduced verticalsshedr. Coup ledawgeeiethemertfects Of ssi ae. 
friction in the lower troposphere, short waves (wavenumbers 
12-15) stopped growing earlier than longer waves which 
extended throughout the troposphere. Gall (1976b) concluded 
that wave intecaceim with the zonal flow plus surface 
friction modified the waves and the zonal flow SO. uilicee 
wavenumbers 5-7 dominated the flow. 

It was observed in the experiments with the genesis jet 
in the previous section, that after G40 had developed, other 
disturbances somewhat larger than wavenumber 40 developed. 
It was suspected that a maximum in the growth rate spectrum 
existed between wavenumber 8 and 40 for a slightly modified 


genesis jet. Therefore, in this Secumen sumer chirectsmon 
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wave-zonal flow interaction are studied for the mean winter 
mmgegenesis jets. These experiments (set B) differ from the 
meevrOus experimencvs (set A) in the following way. First, 
the model is integrated in time until the linear growth rate 
meee Waves Ceases. This occurs at 120 h for the mean 
winter jet and at 84 h for the genesis jet. Second, the 
mitcemeans are computed for u,v, T and from these fore- 
Sareervelds andmare used as the initial fields. The growtn 
of the waves in the modified mean winter and genesis jets is 
then calculated exactly as in the set A experiments. 

ie Growth rates for the set B experiments are compared 
in Fig. 3-15 with those from set A. For the mean winter jet 
modified by the wave-zonal flow interaction the growth rates 
are reduced by 10-30% for wavenumbers Gredteonr =enaneo. | [here 
1s a tendency for the growth rates of the highest wave- 
numbers to nave the greatest reduction. Very significant 
Changes occurred with the modified genesis jet. Here again, 
growth rates are reduced for all wavenumbers greater than 9 
with the most significant reduction at the shortest wave- 
lengths. In fact, the growth rate for wavenumber 40 
decreased from 1.82 day-! to 0.82 day-!. In the previous 
PJeeloneie Was shown that this wave involves both baroclinic 
Mmiseanuiity and Kelvin-Helmholtz instability. Stone (1966) 
and Orlanski (1968) demonstrated that R; must be < 1 for the 
Peer to exist andes tnis Gendition was satisfied (cf, Fig. 


3-2) for the set A experiments. However, tne modified 
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Genesis jet for the set B experiment does not meet this 
eondition. R values in the vicinty of tne modified 
genesis Jet in the lowest levels are in the range of 2-5 
Seer s-lo), which is a small but significant increase 
@emoenve tO the set A experiment. This R. inc redsem)| Ss due 
momame increase of the static stability in the lowest levels. 
The decreased growth rates for the shortest waves are 
consistent with the increased R. values (Stone, 1966; 
Meeeaenski, 1968; Staley and Gail, 1977). R values poleward 
of the jet in the lowest levels are unchanged from the set A 
experiments and remain less than 2. 

The growth rate for wavenumber 20 for tne genesis jet 
meroeetlg Ee redticed as much as the other waves. In fact, the 
growtn rate spectrum for the modified genesis jet exhibits a 
maximum at wavenumber 8 (1.17 day”! ) and also at wavenumber 
0mm 1.2 a). It is suggested that the peak at wavenuniber 
20 did not appear inthe A experiment due to the magnitude 
Gm ene Growtn wate for the “mixed baroclinic mode. 

To study wavenumber 20 (G20), the model was initialized 
Fan@ueeyust this wave in the domain. Here, a horizontal 
mesovucion of 50 km with 28 intervals in the x-direction 
and 128 intervals in the y-direction is used. The initial 
fields were generated as in the B experiments described 


earlier. 
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Fig. 3-16. Initial Zonal ly seveiogeaerneniard Som number 
for modified genesis jel .meommotin sive r.cueees 
2. Values above 10 are not contoured. 
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The time evolution of CE, HCK and VCK for G20 is shown 
eee o-i/. erm CZ 1s dominant, which indicates that 
this wave receives energy through conventional baroclinic 
Baocecees. erm VCK is also positive but is=only one-fourth 
Gf €E-™ Term HCK is negative and about one-half the magni- 
tude of VCK. 

Am @asSt-west cross-section of vertical motion, tempera- 
aieememia trougn position for G20 is given in Fig. 3-16. The 
trough tilts westward with height with the largest slope 
between 850 and 700 mb. The warmest air extends from tne 
mecee trough position to 750 mb with little or no tilt. 
Above 750 mb, the warmest air tilts to the west. With the 
mepongest Vertical motion occurring just east# of the surface 
MeougumpoSsition, generation of KE 1S a maximum in the lowest 
eee ims pdttern is similar to that observed in G40 in 
impeGmerous section and is consistent with the shallow. 
nature of G20 and G40 compared to W12 and G8. 

The surface pressure pattern for G20 is shown in Fig. 
3-19. A well-defined surface low develops about 1000 km 
Glomtn Of the jet stream. Although the horizontal scale of 
eases smaller than that of extratropical, synoptic-scale 
waves, surface pressure patterns in both waves are similar. 
In addition, the energetics for both waves are similar witn 


large positive CE and smaller, negative HCK. 
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Fig. 3-17. Similar to Fig. 327 excememoienanse immerse” 
in modified genesis jet. 
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Surface pressure at 84h for wavenumber 20 

in modified genesis jet. Contour interval is 

4 mb. Figure represents a subsection of the 

forecast domain.  SowUthern@peundary of f ice 
1s at the center of tne forecast domain. 

100 km of right side of domain is the same 

as leftmost 100 km due to periodicity. 
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Striking similarities exist between wavenumber 20 for 
Piemmoaditied genesis jet (G20) and the observations of 
Pacific polar lows by Reed (1979) and Mullen (1979). Reed 
(1979) observed that the polar low has typical wavelengths 
from 1000-1500 Km. Wave G20 nas a wavelength of 1400 km. 
Mullen (1979) presented composited data which exhibited R 
Values of 1-2 at the position of the polar low from the 
Surface to 850 mb. This has also been found to be the case 
Mmemmacod (ct, Fig. 3-16). Mullen (1979) also showed that 
eqempoldr !ow formed 500-700 km north of the jet while G20 
momnecmanbout 1000 km north of the jet. Differences exist 
between the vertical structure of observed polar lows and 
G20. Reed (1979) documented a Pacific polar low which 
ero rced geopotential! perturbDations at 300 mb as large 
as those observed at the surface. G20 has been found to 
exhibit maximum amplitude at the surface with little 
evidence of the wave above 500 mb. It has Deen snown by 
Gall (1976c) that latent heat release can serve to deepen 
~weoreems- iIis process will be examined in Shapeas ee 

It has been demonstrated that G20 is formed by conven- 
Siemaneodroclinic instability which supports the hypotheses 
of Reed (1979) and Mullen ES) that Pacific polar lows are 
Domoewmrecally Unstable disturbances. Energy conversion | 
processes suggest Kelvin-Helmholtz instability is a con- 
aoe Ima factor, DUt 1t 1S dbout Tour times smaller in 


magnitude than CE. In addition, the Kelvin-Helmholtz 


fl 


Instability creates a Structure different Trem thenoosen es 
structure of polar lows. Wave G20 1s BDangtyopiGal ll) saaieee 
5 \WCCrar ee GuUr 

The results for the genesis jet disagree with the 
findings of Sardie and Warner (1983). Using a 3-layer 
quasi-geostrophic model, they found the most unstable wave- 
length to be 2700 km (about wavenumber 10) with a growth 
rate of 1.29 day-l. They did not find large growth rates 
for tne "mixed" baroclinic wave (G40) because quasi- 
geostrophic models neglect the non-geostrophic effects 
necessary for Kelvin-Helmholtz instability. Sardie and 
Warner (1983) showed that inclusion of moist processes 
shifted the most unstable wave to shorter wavelengths (0600- 
2300 km depending on the heating distribution) and concluded 
Chat Men sie baroclinic processes are the primary cause of 
Pac it We mele; ones 

The growth rates computed in this study for the modified 
genesis jet indicate that moist processes may not be neces- 
Sary for polar lows to form. Mullen (1982) documented a 
case of deep continental polar low development with a scale 
of 1800 km in a domain with negligible water content. The 
formation.of a polar | OWwowitin ti 1S sienel 2 cipal eeaiie 
without moist processes 1S .COnS 1 Stent sw ieee findings olnr 
this study. Some other process(es) (possibly terrain 
induced) must be responsible for the deep structure of tne 


continental systems observed by Mullen (1979). 


ie 


DemeeerrecC?l OF FRICTION 

Sieex Per imemtuiset C) Similar to those described in the 
previous two sections is performed to study the effect of 
friction on tne linear growth rates. These C experiments 
differ from the set A experiments in that the surface 
Meer ot 1S imeluded. he zonal mean quantities of u, v, T 
and m are again restored each time step. 

The incorporation of friction in the C experiments 
G@ifiters from the form described in Appendix A. Here, the 


Siienrace stress 15 
—T = oC) IV |v or) 


where Cy = 1.0x10°°. It is assumed that the surface stress 
decreases to zero linearly with pressure through the lowest 
meommodel layers (approximately 120 mb). This formulation 
ere Surfdce friction is used to provide linear friction 
Memetnie entire grid. Although this form is proportional to 
the velocity squared, interactions with the zonal flow are 
zero since the zonal flow is held constant in time. The 
interactions between waves 1S small since the perturbation 
is) scams 

The major effect of the friction is to decrease the 
growth rate of all the waves (Fig. 3-20). Also, the 
Shortest waves are most highly damped, which 1s consistent 
with Gall (1976b) and Gall et al. (1979). This is expected 


Since tne short waves have been shown to be shallow and to 
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nave their largest amplitude in the lowest level. Growth 
rates for wavenumbers 6 and 9 for tne mean winter jet and 
wavenumber 6 for the genesis jet are reduced by less than 
10%. All other growth rates are reduced over 15%. Growth 
rates for wavenumbers greater than 18 for the mean winter 
jet could not be calculated because these perturbations were 
strongly damped. 

Tne most unstable waves for tne mean winter jet in the 
C experiments is wavenumber 9 witn a growth rate of 0.48 
day~!. Therefore, the effect of friction on the mean winter 
jet is to snift the most unstable wave to a longer wave- 
length (wavenumber 12 in set A versus wavenumber 9 in set C 
experiment). For the genesis jet, a maximum occurs at wave 
number 8 (1.04 day-1l) as in the set A experiments and at 
mevenimper 36 (1.38 day-l).§ The growth rate spectrum fot 
the genesis jet is very flat for wavenumbers greater than 
Paeewnhich suggests that the helical rolls formed in the 
“mixed" baroclinic wave tend to be larger when frictional 
effects are included. 

The growth rates for the modified genesis jet with 
linearized friction (experiment D}) are shown in Fig. 3-21. 
poids, the shortest wavelengths are most highly damped with 
the growth rates for wavenumbers greater than 28 too smal! 
memeaiewlate. ihe peak in the growth rate spectrum is 0.99 
day-l at wavenumber 8. This indicates that the combination 


vee sonal tiOw interaction and surface friction alter 
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Gre Growth rate Spectrum of the original genesis jet such 
Cnat tne synoptic-scale wave dominates. Clearly, tne deve- 
lopment of a polar low must involve processes other than 
moet inic instability to attain observed intensities. The 
effect of diabatic processes on the polar low will be 


studied in Chapter IV. 


E. CHAPTER SUMMARY 

The experiments described in this chapter lead to two 
mmpOrtant points. First, significant differences exist in 
Miemd@owlm rate spectrum for the two jet profiles. It was 
found that the mean winter jet yields a maximum growth rate 
fOr wavenumber 12 (0.59 day-!). This result supports the 
findings of Gall (1976a), who computed linear growth rates 
fOr a mean zonal flow similar to the mean winter jet. He 
found wavenumbers 12-15 to be most unstable with a growth 
Puen 0,0 day-l. [n addition, the shape of the curve in 
migeec-| tor the méan winter jet is quite similar to that 
given by Gall (1976a), who included only wavenumbers < 15. 
The genesis jet case yields a maximum growth rate for wave- 
mambicnmeoetl.,22 day-!) and 40 (1:82 day-1). It was shown 
tnat wavenumber 12 for the mean winter case and wavenumber 
8 for tne genesis case are baroclinic waves which extend 
throughout the troposphere. However, wavenumber 40 was 
found to be a baroclinic/Kelvin-Helmholtz wave which 1s 
confined to the surface-/50 mb layer and is characterized 


Deeceiietica! type circulation about a Morizontal axis. 


a 


The second major point is that the linear growth rate 
Spectrum 18 strongly attiected by process com wind Gemmce 
continuously modify tne initial zonal flow. The processes 
tested here include wave-zonal flow interation and surface 
friction. I[t was foumd that SUmieice si tice el cao ceece 
wavelengths, Dut is especially effective in damping the 
shortest waves. In daddijriron..sit,,a¢eu rp ereemon shifted tne 
most unstable wavenumber for the mean winter jet from wave- 
number 12 to 9. It was also found that wave-zonal flow 
interaction modified the shortest waves the most and led to 
a significant change in the Growin rawemopec ry tpt Oo mmeme 
genesis jet. The baroclinic/Kelvin-Helmholtz wave was 
damped considerably, as the growth rate was reduced froin 
1.82 day-! to 0.82 day-!. This was due to an increase in 
the low-level R values from less cham to Va lilies Melee =e 
The modified genesis jet exhibited maximum growth rates at 
wavenumbers 8 and 20. It was found that wavenumber 20 was a 
baroclinic wave with maximum amplitude in the lowest level 
and resembles tne Pacific polar low described by Reed (1979) 
and Mullen (1979). In addition, this wave is not as deep as 
the longer baroclinic waves, with the maximum upward motion 
occurring at 800° mb as “Opposed tows) We fon cite onde: 
waves. Friction 1s found to damp wavenumber 20 strongly, 
which suggests that polar lows cannot attain observed | 


Intensities throtigh baroeynyceams Gen ey ecm eotes 


iS 


IV. NONLINEAR MODEL SIMULATIONS 


~- 


It was shown in Chapter III that the linear growth rates 
for both the mean winter and genesis jets are largest for 
relatively large (wavenumber 8-9) baroclinic waves when 
Peretional effects and wave-zonal flow interaction are 
considered. In addition, tne genesis jet exhibits large 
peeve rdtes for two smaller-scale systems. The first, 
wavenumber 40, was found to be very shallow and forms 
mipougieparoclini¢c and Kelvin-Helmholtz instabilities. After 
this wave modifies ene zonal flow in the vicinity of the 
jet, a peak in the growth rate spectrum occurs at wavenumber 
Zee ee km). This particular wave, G20, forms through 
conventional baroclinic gaa lacy and resembles the 
Pacific polar low described by Reed (1979) and Mullen 
(1979). However, the vertical scale of G20 is smaller than 
observed for the Pacific polar low. 

The development and structure of G20 will now be studied 
in detail. This will be accomplished through the examina- 
tion of model output and energy budgets of tne storm at 
various stages of development for both adiabatic and 
diabatic simulations. The model configuration to study 
mmesewwaves is the same as that used in studying the energy 
transfer in these waves in Chapter [II {Table 3-1). A 


similar study of wavenumber 8 for the genesis jet 1s 


te) 


presented in Appendix C. The results of the experiments in 
Appendix C are used Tom Compan i someo) Mewe pola, 5 lo 7secascr 
CXC raiamom ec | sy 10 oti CMeCene system. 

The following procedure is used to incorporate the 
initial perturbation into the fields for the model simula- 
tions. First, the perturdetion@™ for 4a single Wave 1s impo 
in the initialized fields in the manner described in Chapter 
II (i.e., constant amplitude in the ronan ft Second, the 
model is integrated for 36 h with the initial zonal mean 
values of u, Vv, 1, 9 and wt restored every Cimemovepeunoce 
inviscid, adiabatic Comd?t1o0ms. 1 11sec to W Suma nomena 
reach the point at which exponential growth starts as in tne 
Chapter [If experiments. This insures the proper vertical 
structure of the WaVe at the infttaecqme. “(ie gem oe 
forecast fields serve as the inttral fields ror alle simncm 
tions presented here. 

The sea-surface temperature for these experiments is 
computed as follows. The mean air-sea temperature differ- 
ence given by Mullen (1979) for the genesis conditions is 
interpolated to the forecast grid in the same manner as for 
tne u-component described in Chapter II. Mullen (1979) 
found the maximum air-sea temperature difference of 6.8°C to 
be approximately 750 km north of the jet axis and decreased 
to about 4.5°C to the north and vodpod wee sees Ommrimmen ernie 


jet axis. This difference is added to the lowest level 


initial temperature field. The sea-surface temperature and 


30 


air-sea temperature difference fields are shown witn the 
individual experiments. Surface vapor pressures are the 
Saturation values corresponding to the sea-surface tempera- 
tures. Sea-surface temperature and moisture values remain 
PomSeant over each forecast. 

The diabatic runs use the complete model while the 
adiabatic runs neglect precipitation physics and sensible 
Simaemreacent Meatingmwithin tne PBL. Also, the adiabatic runs 
Peaaei4eerne Siitmace momentum fluxes to incorp@rate surface ~ 
Miieronmeane the surface heat flux to predict the PBL deoth 
even though the sensible neating due to the surface heat 


flux is not added to the temperature field. 


A. WAVE GROWTH AND STRUCTURE 

pone on Or ume minimum central pressure versus time of 
620 for both the adiabatic (G20A) and diabatic (G20D) simu- 
ememens 15 given@im Fig. 4-1. The inclusion of diabatic 
Cffects results in G20D attaining a central pressure of 
Poe Whicm 1S 19 mb lower than GZOA. Clearly, diabatic 
Seeieese serve towintensity the disturbance. This 7s 
expected based on the studies of Bullock and Johnson (1971), 
Tracton (1973), Gall (1976c) and Sanders and Gyakum (1980). 

Mullen (1983) computed deepening rates for several 
Meetric 0Glar lows accoOrdimgs to the method? proposed by 
Sanders and Gyakum (1980), who defined 1 Bergeron (8) as a 


Bemeealepressure fall of 1 mb h-! for 24 h at 60° latitude. 
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Ine deepening rate at any latitude 4, is adjusted by multi- 
plying by sin 60°/sin 9. Here G20D exhibits a 13 mb fall 
Meo mean o2Ne which G@orresponds to a deepening rate of 
ioeen es ine maximum deepening rate for G20A is 0.22 B. 
Mullen (1983) showed that Pacific polar lows can achieve 
@eepening rates of about 1.9 B, although this seems to) 6ccur 
only when the polar low interacts with the main frontal band 
equatorward of the deepening low. For a polar low which 
did not interact with tne main polar front, Mullen found a 
Geepeming rate of 0.85 8B. This corresponds well with the 
deepening rate computed here for G20D. Observed deepening 
rates, such as those computed by Mullen, should be expected 
Pomeeceidarger than those computed in this study. Pacific 
polar lows generally deepen as they move into a large area 
of low pressure often centered in the Gulf of Alaska. 
Therefore, some of the pressure fall must be attributed 
HOMmiomnew environment of the storm, rather than actual 
miecwsii ication of the individual cyclone. This “change of 
environment’ did not exist in G20D and therefore, deepening 
rates computed using real data could very well be larger 
tnan tnose computed here. 

The sea-level pressure fields at 60, /2, 84 and 96 h for 
G20A and G20D are shown in Figs. 4-2 and 4-3, respectively. 
During intensification, the low center in G20D (hereafter 
G200) moves due east while G20A moves toward the east- 


Momeredse. At 965 Nh, the center in GZ0A 15 about 200 km 
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represent a subsection of forecast domain. 
Southern boundary of figure is located 200 km 
northgof center sof Tetgeeast domain. See kines 
right side of domain is the same as the leftmost 
100 km due to periodicity. 


84 


15 


OISTANCE (101/108) 


10 15 5 


DISTANCE (HH/100) 


ss) 


Peds. 


3 ld 
OISTANCE (KM/12@) 


4-3. 








Sout) aie eo les 


2) 0) 


15 


ee 


OISTANCE (iV! 


RY 





15S 


DISTANCE (KIN/1001 1 


5 





2) 18 
DISTANCE (KM/12@) 


Game 4=c2 except for G2Z0D. 


north of the corresponding position oT G2ZUD 2) Ss iice a ice eee 
processes are included in G20D, it is Suspected tnat the 
sea-surface temperature (Fig. 4-4) plays a role in deter- 
mining the position of the low. Sanders and Gyakum (1980) 
demonstrated tnat strong oceanic development tends to occur 
near the maximum gradient of the sea-surface temperature. 
Here, it is clear that during the developing stages, the 
center in G200 remains in the region of maximum gradient of 
Sea-surface temperature (near the 6°C isotherm), while the 
center in G20A moves away from it. 

The low in G20A is observed to form about 900-1000 km 
north of the jet axis While tne low 1p Gator orice 
750 km north of the jetmrekIS., [his 1S GonsSistenmt wi caeee 
composite results of Mullen (1979), which indicated that 
the polar low forms =ogGg—=7ae fac Porth of tne jet. Therefore, 
realistic deepening rates and more realistic storm DigsS 1ew Cline 
are produced when diabatic effects are included. 

Both G20A and G20D) exhibit a Band of positive wearne- 
lated u’ and w Similar to that Gin G4Q)( cf er oc eem 
parallel to the trough extending from the low center. In 
this region, the R- < 1) Wii Gnis Ug giesic seed one aie 
Helmholtz instabdility May be S19mit tea mem mie die cuee olin 

East-west cross-sections of helgnu perGurode lomse atc 
are given in Fig. 4-5 for G20A and G20D. JInese figures are 


representative of the developinu stages of the disturbance. 
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G20A exhibits a maximum deviation at the surface of about 
Pewee (ae Wave is damped strongly in the vertical so that 
maieeampiittude at SOQ mb is only 50% of the surface value. 
G20D is noticeably stronger with a surface maximum amplitude 
of 110 m. The amplitude does not decrease as rapidly in the 
Poeeted! 2S in G20A" The amplitude at 500 mb is about fone- 
nalf the surface amplitude. 

smie@e G20) extends throtigh a greater depth than G20A, it 
should be expected to have a phase speed greater tnan that 
of G20A. Phase speeds computed from the surface pressure 
eemeencein Figs. 4-2 and 4-3 are 12.4 m s-! for G20D and 
Pome! tor G20A. Somleman the increase in @hase speed of 
eeubemust a41so be due to it being 1mbedded in stronger flow 
weomocne more southern track. 

Best -west cross-sections of Vertical motion and temper - 
ature at 84 h are shown in Fig. 4-6. The reason for the 
increased vertical scale in G20D becomes obvious upon 
iempeerron Of thew» dnd | fields. For GZ0A, the warmest 
Seemex tends upward from the surface to about 750 mb with 
Meee or m0 tilt. the Maximum vertical motion occurs just 
Sasa r wie Warmest dirk which favors generation of KE. 
Sew oeenmbD, the Warmest dir tilts to the west with height 
away from the maximum vertical velocity similar to the 
pattern exhibited by wavenumber 20 in the linear experiments 
in Chapter III. In addition, the temperature perturbation 


meepeases above 750 mb. these factors are consistent with a 
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Gomieca Generation or even a destruction of KE above 750 mb. 
imesew iil! Oe discussed fUrtner in the next Section. The 
aeeucrine Of GZUA Stronglyunesembles that of G20 in the 
meecmme periments im Chapter [1l. This 15 not surprising 
Since no other waves are included in these studies, and 
Since the disturbance is relatively weak, wave-zonal flow 
Mieeecct ion 1S Sitall. in contrast swo the vertical structure 
of G20A, the low in G20D exhibits a strong positive tempera- 
ture maximum which extends from the surface to 450 mb witn 
eine Or NO tajlt. this results ain a generation Of ‘Keewove tr 
a mucn deeper layer than in G20A. 

One reason for the increased depth of the temperature 
perturbation in G20D can be seen by examining Fig. 4-7, 
which shows the amount of total , convective and stable 
precipitation Pndtwaecumulatied diring the 12 h period ending 
at 84 nh. Mr ated Of Precip BLatmon (about 22 mm 15 found 
at and just ahead of the path of tne low center with a 
secondary maximum (about 15 mm) south of the low along the 
PeaG@emm ADOUL /02 of thesprecipitation accurring near the 
Mewmcenter and aqwoeur 854 in the Southern@imaximum is convec- 
dee. Based on cross-sections of 8, and relative humidity 
(not shown), conditional instability extended up to about 
Temibewiile stratiform clouds (relative humidity > 100%) 
Poasenideasie to about 350 mo 1m the warm sector ahead of the 
low. The release of latent heat over this depth is consist- 


ent with the warming observed in G20D anead of the low. The 
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DZ 


diabatic heating effectively reduces the amount of adiabatic 
Seeing In the region Of ascent. 

To separate the effects of latent heat release by 
Peeeipicdtion trom the PBE diabatic processes, a Simulation 
(G20E) was made which was similar to G200 except that latent 
heat release Dy precipitation. was NeewaiNwowed. (GZ0E attains 
minimum Central @pressure of 993 mb (Fig. 4-8), whichis 
moeme lower tha tme adiabatic simmlation, but 9 mb higher 
mame tie tUllyM@dyabatic simulation. Obviously, latent heat 
@oease DY reel pieation 1S note actaing alone in tne intensi- 
feedeyon OfeG2OD. BSignificant contBibutions to the 
intensification of the polar low are made by three otner 
meocesseS. First, the sensible heating due to the large 
air-sea temperature difference insures a strong meridional 
Now=level temperature gradient (i. e., badroclinity) over a 
longer time period than found in tne adiabatic simulation. 
Second, the sensible heating within the PBL can be realized 
in tne mid-levels almost immediately through dry convective 
gagustment. This 1S because the lapse rate is nearly dan 
adiabatic in the genesis jet from the surface to near 750 mb 
iene ev icin ity Ofmmne spo0lar sii. Third, athe large, initial 
ajr-sea temperature differencegacts to maintain a nearly dry 
adiabatic lapse rate in the lowest troposphere (i.e., low 
R.) whicn should increase the growth rates of the snorter 
waves (Staley and Gall, 1977). To show the influence of tne 


surface heat flux, G20E was rerun without a perturbation in 
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tne initial fields. The top of the low level dry adiabatic 
layer near the latitude of the polar low rose from near 

Pee owimicidl ly toeneanso50 5° mb by 96 ume The effect of 
these two processes on the energetics of G20D will be 
aeseussed further im the following section. 

Piewci rece Olmelorardbabrc OroOcemsces on G20 can be seen 
in Figs. 4-9 and 4-10, which show the heignt and temperature 
merase at Of fh TOmecogn 700 and 500 mb fT@r G2Z0A and G200, 
Mespectively. the amplitude of G20D is stronger at all 
Mevels. in addition, the Wanmest (coldest) air in G20A 
remains just to the east of the trough (ridge) position at 
all levels. Jhe amplitude of the thermal ridge jis larger 
Gumoerartmer to the east of the trough axis in G20D, espe- 
exeuy at (00 and S00 mb. 

nice stemace  SensiDle dndy watent heat fluxes for G20D at 
EOoumeare shown if Fig. 4-11. Tne maximum for both is about 
Somiem-csand 15 found to tme southwest of the low eet 
Deviate cOlad front. This, Occurs because this quadrant 
exhibits relatively low temperatures and the strongest low- 
level winds. Although downward fluxes are typically 
observed in the warm sector ahead of alow, G20D exhibits a 
secondary maximum of upward surface fluxes roughly one-nalf 
the fe qa values in the cold air. There are two reasons 
Toa emp widrd Ghlixes ein tie Warmsdir. First, the initial 
air-sea temoerature difference in G20D was at least 5°C in 


the vicinity of the developing low, a condition that results 
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in an upward surface heat flux in all quadrants around the 
Tow. Second, the low-level norizontal temperature perturda- 
tion associated witn G20D (Fig. 4-10a) did not become strong 
enough (at least through 84 h) to reverse the sign of the 


initial air-sea temperature difference. 


Pee ENERGETICS 

Tne equations of the energy budyets used in this section 
are given in Appendix B. However, the terms involving 
mearmmoges In tne stirface pressure (1. e., PAE, PAZ, PKE, PKZ) 
are small relative to the other terms, and therefore are not 
Miemtaed in this discussion. Ihe budgets are not perfectly 
poememeed die tO Eruncation, interpolation, round-off errors, 
eeueor the adiadDatic runs, the residual is generally an 
order of magnitude less than the leading term(s), and the 
Pood) FOr the diabatic runs 1S about 5 times less than 
the leading term(s). 

KE Dudgets, integrated over the forecast domain, are 
presented in 6 h intervals in Tables 4-1 and 4-2 for G20A 
and G20D, respectively. Steady increases of KE are noted in 
munis 7rom the beginning of the forecast. Tne maximum 
ameceraned in G20A (7 x 10S J m4) 1s an order of magnitude 
eominan the maximum reached in G20D (8 x 10° J m-¢). In 
addition, the maximum value occurred from 114-132 h for G20A 
mies culn fOr G200, which shows the importance of tne 
Maleate for ecesses im increasing Doth tne deepening rate as 


Pelseedsetne Intensity Of the system. 
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Time-heighnt cross-sections Of KE TOr GZURRandeG le ae 
given in Fig. 4-12. Por G20A, the maxi mumere soce ur sea 
mb and decreases rapidly WaTOnen Git: For G2Z0D, the marin 
KE is found at 90 h, also near 950 mb, but a secondary maxi- 
mum of about 2/3 of the 950 mb value exists near 450 mb. 
Tnis 18 a reflection of the ineweased deweio i e011 seus 
as moted in the previousmseccule TT 

During the growth of G20A and G20D, CE is the largest 
term, which indicates that baroclinic conversion is the 
principal mechanism in forming the 1400 km wave. The values 
of CE are more than an order of Magnitude larder fom Glu 
than for G20A. Time-neight cross-sections of CE are shown 
in Fig. 4-13. For G20Ay CE hasta maximum ct el20n enn 
875 mb with negative values above 700 mb as tne storm 15 
developing. TnNtS 15 GOns1Stene Mya amie men le cl; ec men 
G20A found in the previous section. It was noted tnat a 
thermally direct circulation was eVident below about / 500i 
which 18 consistent with the shallow atte OfeCr 1 aeiage 
4-13. For G2005 CE exhibits 4 maximum abt e4 fnear 7 Ocean 
that 1s nearly an order of magnitude larder than the G2e 
maximum. in addition, positive Valwes so te eee. bona. 
nearly 400 mb. It was noted in the previous section that 
the effect of diabatic heating was to shift the warmest air 
above 750 mb to the region of the strongest i pwan cee cee 
thereby enabling the generation of KE, through the term CE, 


to strengthen the circulation in the middle to upper levels. 
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Tne norizontal shear, or the barotropic term, HCK, is 
about one order of magnitude less than CE for both G20A and 
G20D during the early stages of development. As tne wave 
grows, HCK is negative and increases in magnitude. In G20A, 
HCK reaches a value of about one-half that of CE at 120 h 
marie in G20D, lUkeibecomes tne largest term at 96 h. Inere- 
Meme OdrOmnopie effects are important in halting the growth 
of the polar low, Ssoeecnaliem when dlabdatuterwet fects are 
included. Time-height cross-sections of HCK are given in 
Fig. 4-14. As with CE, G20A exhibits a maximum of HCK near 
Pema wit le G20Dse ani bitsee maximum at 850 mb and near 
SO Gein . 

Tne Kelvin-Helmholtz term, VCK, is positive at all times 
gov BOLI TOrecasts. It was snown in Chapter rll Unaoaeee 
involves the correlation between u' and w' multiplied by the 
vertical wind shear. During the development of the polar 
low in these experiments, strong positive correlations were 
observed between u' and w'. Generally, this implies 
easterly, upward motion ahead of tne trough and westerly, 
downward motion behind the trough. Coupled with the large 
Mercicad| Wind Snear, large, positive values of VCK result. 
It was shown in the previous section that the sensible heat 
flux in G20D reduces the stability in the lowest levels, 
pimemeresultsS mama decrease Of Ri> Which is consistent with 
the large values of VCK. However, VCK is ruled out as the 


term which is responsible for the rapid deepening of the 
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Poo eow. in fact, VCK May inhibit polar low formation. 
As described in Chapter III, positive values of VCK lead to 
Memrealscirewlations mn thesvertical plane whicn are quite 
Meirerenc trom the Observed Structure of polar Lows. The 
presence of this helical circulation will tend to decrease 
Siem vertical wind shear (7.6., baroclinity) and eit in 
mescaticed Growth rate of mame polar low through CE. Also, 
ing the rapid Growth of G200 from 78-90 h, GE is the 
dominant term. For example, VCK is threa times less tnan 
CE at 78 n. When weaker development was occurring prior 
momo VCK Was banger titan CE ingGZ20D. 9 In G20A, VCK was 
teiemally ad faetor of two smaller than CE at all times. 
Mme=melgnt cross-sections of VCK are snown in Fig. 4-15. 
Poa, the increased depth of the disturbance is reflected 
by a double maximum in VCK at 850 mb and 500 mb in G20D and 
only a low-level maximum of VCK near 850 mb for G20A. 

ice dwar Oleol KE DY friction, DE; 1s negative at 
eeebimes, wien indicates a sink of KE. For G20A, DE 
slowly increases in magnitude during the forecast so that 
meee, Dewseone half that of CE. At 114 h, DE is larger 
Mimo wich snows the strong Influence of friction in 
Sao unvoeeem)S disturbance at all times. in fact,.DE is the 
largest term in the G20A budget from 114-138 h. This is 
expected due to the shallow nature of this disturbance. 


Be wGcUUN DE Dlays ad reduicedepuemscil mesyqnificant role. 
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Peenoegneve atiains d maximum at 90 h in G20D when KE is a 
iaxwMum, DE 1s still 30% smaller than HCK. Therefore, when 
diabatic processes are included, barotropic damping is 
Sominant 1n halting wave growth. This is due to the 
Mitem~eased depth of G2Z0D. Frictional effects are relatively 
small in G20D since they are felt only in the lowest levels 
of the model. However, it was shown earlier (Fig. 4-14) 
Midibeion exhibits ldrge values near 850 mb and near 500 mb 
ana 1S, therefore, more efficient in reducing overall wave 
Growth. 

AE budgets, integrated over the forecast domain, are 
presented in 6 h intervals in Tables 4-3 and 4-4 for G20A 
guemacuu, respectively. G20D exhibits more AE than G290A 
at all times. The maximum AE for G20A (2.6 x 103 J m-2) 
imeround at 114 h, and Eee Maramum AE tor G2Z0D (1.3 x 
oe J To) VSmepolnd eat oO) Hh. meiherefore, as avidanaee with 
KE, AE develops faster and reaches larger values with the 
Miemesion Of Gdlabdtic processes. fime-heignt cross-sections 
Of AE are given in Fig. 4-16. The major difference between 
mimemporecastsy1s thdt G20D exhitbD1ts a strong maximum of AE 
emo) MD. §)N1S Teadture iS™a0sSent In the G20A forecast. 

The generation of AE by diabatic processes, GE, 1s 
positive at all times in G20D asa wom ee ee 
monet ound Tm davance of the storm center (cf, Fig. 4-8). 
Mtemedwidescavdnues Of GE occur Trom /S8-90 h during the rapid 


mimenmsitedtdomeor G2OD. The time-heignt cross-section of 
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meetor G20) 1S shown im Fig. 4-17. Maximum values of GE are 
found near 800 mb with positive values extending from near 
950 mb to 400 mb. Small negative values are observed near 
mnemoimerace Gwe tO the largest surface fluxes occurring in 
eiemwaoid air behind the low (Fig. 4-11). 

lime norizonta!l eddy Wheat transport term, HCA, 1s 
positive during the development of G20A and G20D, whicn 
indicates a generation of AE. Maximum values of HCA for 
G20A are found near 950 mb with smaller, negative values 
cere U0 mb (Fag. 4-18). For the diabatic run, HCA is 
larger in magnitude and the positive values extend up to 
yom. 81nis tllustrates thaw the diabatic processes have 
Pt eanin an incweased@errem dation above 750 mb such thet 
norizontal eddy heat transport can generate significant AE 
at these levels. This is expected due to the height and 
temperature patterns exhibited by G20A and G20D (Figs. 4-9 
and 4-10). 

micewverticawecomponemy of CA, VCA, 1s also posi@ive 
during the development of the polar low in both experiments. 
SaamieeGacOD EXMIDIts ValUes Over one order of magnitude 


PerdearmitiiatmacUA. es | Chougm the term congaining ee favors 





Ceealict ion Or AE, UmemtcrmSmeomt aining w E' 3f/ 9p and 





w' T'3E/3p are positive and larger in magnitude, and result 
in an overall gain in AE through VCA. Time-height cross- 


sections of VCA for G20A and G20D are given in Fig. 4-19. 
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iiemmest Stemificant differences between the forecasts 1S 
mie dange weredase 1n VOCA near 600 mb in G20D, whicn again 
reflects the increased circulation in the middie layers. 

The results of these energy budgets illustrate now the 
aed tic effects ledd to an tncrease in storm intensity. 
tiene tne rapid intensification of G20D from 73-90 h, the 
wi eogencAwand VCA \isSe6-=10 times larger than GE. [nis 
Suggests tnat mucn more AE is formed by eddy heat transport 
processes than generated directly by diabatic effects. 
ims 1S supported by the fact that Sarin G20, w hase 
neglected latent heat release by precipitation resulted in 
a1l10mb decrease in minimum central pressure from the adia- 
weeersimuldtion. inerefore, 16 1S concluded that diabatic 
processes in the Pacific pe lame 0 riiiar i ibeeserve tomd|lter 
tne overall temperature structure of the disturbance such 
met eecdgdy neat transport, and therefore, AE values are 
[mmemedsed. Ine rapid storm growth is due to baroclinic 
Som vension Of this AE to KE. 

The results presented here agree with the speculation of 
Sapavesand Warner (1983) that moist baroclinity is the cause 
meee POld? |OWSeme ine linear studies in Chapter 111 
‘demonstrated that wavenumber 20 exhibits growth rates as 
large as wavenumber 8 after a period of modification to tne 
Genesis jet. Although wavenumber 20 forms in the nonlinear, 
Tce cmieGeimw1th Trverion, tne intensity of the polar 


low is much weaker than observed. Furthermore, the wave 


tends to be too shallow. Incision Of (G@iabat eer eo 
leads to a significamtemamerease TA Fnite nsec uulecumo 
lower central pressure) and a deeper, more realistic verti- 


Cdili7Struc uure. 


C. EFFECT OF ATR-SEA TEMPERATUR GO rir ee 

In the previous two sections, large differences were 
observed between the adiabatic and diabatic experiments. 
The air-sea temperature differences given by Mullen (1979) 
for the genesis case (cf, Fvg@ieet=4) Were Ue cd rious tiem 
fields. The results of twosaddi?rcional s tmineae rons —camemeree 
are presented here to illustrate the effect of tne air-sea 
temperature difference. Details @f di!l four G20 Ss) nu ioe on 
are given in Table 4-5. A zero initial air-sea temperature 
difference is specified in G20B, white in G20C the arrases 
temperature difference is one-half of that used in G200. 
All diabatic processes used in G20D dre also used in G20B 
arivaesG 2 0ne 
Table 4-5. Physical processes and initial air-sea tempera- 

ture differences for polar lowes Tmiiiac 1 ouee 


Initial Air-Sea 


Diabatic Temperature 
Experiment Frictron Processes Difference (K} 
G20A Yes No Mullen (1979) 
Genesis Case 
G20B Yes Yes 0°20 
G20C Yes Yes One-half value 


used in G20A 


G20D Yes Yes Same as in G20A 
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The minimum central pressures versus time for all four 
G20 simulations are given in Fig. 4-20. The least intense 
storm was produced by the adiabatic simulation, G20A, which 
attained a minimum pressure of 1003 mb. Although tnis 
Forecast did not include diabatic effects explicitly, the 
initial air-sea temperature difference did play a role in 
mime yO|ULTON Otevme StOrm since the surface friction 1s a 
miecrem Ot the Stapility of the PBL (Eqs. A-10 through 
Peeole Initially, the PBL was unstable due to the large 
air-sea temperature difference imposed and remained unstable 
Mimeugnout tne forecast. This is somewhat unrealistic 
because the surface wane PIM mou ndmac tet ocdr ive the P Bat 
Coward a neutral stability. Therefore, G20A was rerun using 
an initial air-sea temperature A eenonence of zero to deter- 
mine the effect of PBL stabilization on G20A. Although no 
difference in storm structure or Speed was noted, the 
minimum central pressure did decrease 2 mb more than in 
Gz OA. 

Pie simulate nemeoef GZUD produces the most intense storm 
and the storm with the greatest vertical depth. Tne rapid 
deepening of the system is crucially tied to the latent heat 
release. From 72-84 h, when G20C and G20D exnibit the most 
Significant deepening, there is 15-22 min Of precipitation 
near and ahead of the storm center in G20D. During the 
Same time, G20C has 5-7 mm and G20B has less than 3 mm of 
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ie) 


ie oreresmlts OF experiments G20B, GZOC and G20D 
mmistrdte the type of environment tnat is favorable for 
Significant polar low development. As the polar low 
develops, relatively cold air iS advected to the west of 
the storm in the lower levels, while warmer air is advected 
ahead of tne storm. This temperature advection drives the 
PomeboOwdrd more unstable conditions behind aie storm and 
Mowdra more StadDle Conditions ahead of the storm. For an 
Piedad! ly netitral PBL {as in G20B), this implies an unstable 
moe tme €Old air and a stable PBL in the warm air during 
mien agrOwtn stage. Iherefore, surface fluxes in the warm 
meebon rend tO be relatively small. This implies that 
momatively little or no additional moisture will be added to 
the system which may be used to sustain additional growth. 
mumeoaalt1On. surface friction would also be relatively smal] 
in the warm sector which implies smaller vertical motion and 
smaller moisture convergence, and tnerefore a reduced amount 
Peeerecioitation. in the unstable PBL behind the surface 
ewes reldvively large surface fluxes exist. However, 
meating in tnis area will destroy AE. More favorable condi- 
Mmeonseroredevelopment exist fOr an initially unstable PBL in 
all quadrants around tne developing disturbance (G20D). In 
merscecase, lard] amounts of heat and moisture are added to 
tne near-surface layers, which is a favorable condition for 


Some seeione ino agadltion, surface friction is increased, 


aa 


which suggests stronger vertical emot 1 Onmandmr ome One 
Increase in precipitation. Although AS Sis ide sie aa 
inthe cold air Dehind the lowe seh is sisesc tee co mipiceiiecrern 


to the generation of AE in the warm sector 
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VY. CONCLUSIONS 


imemmunoose Of this study is to determine the reasons 
for the growth and maintenance of Pacific polar lows. In 
Meanereuldr, the major hypothesis to be tested is that tne 
Poem be polar low forms through baroclinic instability. To 
Pecomomish chis, a S3=dimensional, primitive equation model, 
including convective and planetary boundary layer parameter- 
meaqeons 1S used with two sets of initial data. fhe first 
mee loma jet profile typical of mean winter conditions and 
miemsecond 15.0 Jet profile typical of conditions present 
during the genesis period of Pacific polar lows (Mullen, 
1979). A linearized version of the model is used for 
aecedelc and diabatic simulations of the Pacific polar low. 

The linear studies indicate tnat medium-scale 

(wavenumber 12) waves are most unstable for mean winter 
conditions. However, wave-zonal flow interaction and surface 
friction are found to preferentially damp tne shortest waves 
and thereby snift the most unstable mode to wavenumber 9. 
Energy budgets verify that this wave forms through bDarocl- 
iMimenwinstabllity. These results are similar to those found 
Oy Gall (1976a), who computed the growth rate spectrum for 
wavenumbers less than 15 for mean winter conditions. 
Although this study includes wavenumbers 15-30, growth rates 
of these waves were found to be much less than for the 


longer waves. 


The genesis jet exhibits @ peak in tnewenOy ieee 
Spectrum at wavenumber 3, which also grows primarily tnrougnh 
baroclinic instability. However, the largest grow Giese 
found at wavenumber 40. This disturbance 1s confined to the 
lowest 300 mb of the troposphere and exhibits a helical 
Clreula tion. An energy budget verifies that thnis wave 
grows through baroclinic instability and Kelvin-Helmholtz 
(vertical shear) instabilitywemstome (1966) Vancmi ms eanaran 
(1968) showed that Kelvin-Helmholtz instability can De 
Significant for Richardson numbers less than 2s a coma eae 
initially satisfied near the surface in the genesis jet 
case. 

Wavenumber 40 is found to be sensitive to the initial 
conditions and physical processes. ]nis wave Cxnmoltoea 
relatively short period of linear growth, approximately one 
day, as opposed to several days for longer hae Le es 
Suggested tnat tnis 18 due to the stabilization of "thee lowe 
relative to wavenumber 40, by either wavenumber 40 or a 
combination of waves. For this reason, a modified genesis 
jet is created by allowing perturbations to interact with 
the zonal flow for 84 h, and then computing tne new zonal 
flow by averaging all 84 h forecast fields over latitudinal 
circles. In this modified genes i sejet en eiiemr ne manecison 
number in tne low-levels near the jet axis ranges from 2-95 
and the growth rate of wavenumber 40 (650.62 day > > enim 


is in contrast to the oridinal genesis je pees i Cee oeieer sine 
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Richardson numbers less tnan 1 in the low-levels underneatn 
Meemyec dxXis and d Growtn rdate for wavenumber 40 of 1.82 
ey. «6cdin addition, Marcce sit econ 1s 1 OUmdmnoe rest) t 1N 
a further reduction in the growth rate of wavenumber 40. 
However, nonlinear simulations of a longer wave (wavenumber 
20) indicates the presence of Kelvin-Helmholtz instability 
DyeeriemetormatiOn Of bands Similar to that found in wave- 
number 40 in tne linear experiments. These bands are 
located along tne front and are correlated with Richardson 
numbers less than 1 in the low levels. These low Richardson 
numbers ae a result of a large surface sensible heat flux 
forcing near-dry adiabatic lapse rates in the near-surface 
layers. This suggests that although wave-zonal flow inter- 
Beton dnd Surfdce friction act to damp severely the growtn 
of the baroclinic/Kelvin-Helmnoltz wave, physical processes 
Can create favorable conditions for its existence. 

Growth rates for the modified genesis jet are signifi- 
cantly different from the original genesis jet with one 
growth rate maximum at wavenumber 8 and another of equal 
Size at wavenumber 20. The latter disturbance resembles the 
Pacific polar low described by Reed (1979) and Mullen (1979) 
iment forms apouc LO0G™= km morth of the jet axis witn a 
wavelength of 1400 km. However, wavenumber 20 has maximum 
amplitude at the surface and extends only through about the 
lceaeecee mo Of the trovosonere. This 1S tmconsistent with 


observed mature Pacific polar lows which can exhioit heignt 


lea 


perturbations at 300 mb as large as at the surface (Reed, 
1979p Energy budgets verify that wavenumber 20 forins 
through baroclinic instability. The magnitude of the growth 
rates found in this study for wavenumbers 20 and 40 demon- 
strate that systems With horizontal sealieceo @epieer eae 
polar low (and smaller) are unstable without the inclusion 
of diabatic processes. 

The results of tnese linear studies using the genesis 
jet differ from those of Sardie and Warner (1983) (hereafter 
SW), who showed that the maximum growth rate for the genesis 
Jet is fee? day! at wavenumber 10. SW found no evidence of 
the larger growth rates for shorter wavelengths as found in 
this study. Two reasons are cited for Eis. Furst, Sigua 
a quasi-geostrophic model which neglects the non-geostrophic 
term Cro es) that was found to be™signitreane in tie 
development of wavenumber 40. Second, the model used by SW 
contained only three layers, with the lowest layer being 200 
mb thick, whereas the model used here contains three layers 
within the lowest 200 mb. Tm gehisist udyeretor sl iiear seed 
batic, inviscid conditions, wavenumbers 20 and 40 are 
confined to the lowest 250 mb, Whi chietig@es Usmetiicememic 
vertical resolution use@@™ by "SW 7 Siecle yet ame Ommen oem ome 
these features. 

Adiabatic and diabatic nonlinear simulations are 
performed for wavenumber 20 of the modified genesis jet. 


Although tne adiabatic system intensifies, it remains a 


Saailrow GdisturdDance as In the linear experiments and the 
@eerpeming rdlesis much smallér than typically observed for 
Poem re OOldrelows. Ine Inclusion of diabatic effects 
dramatically increases both the disturbance depth and 
deepening rate. This is accomplished by enhancing the 
temperature eaten tis (eit une dt StUTDance: such tiat 
large amounts of AE are generated by eddy heat transport. 
Tne temperature perturbation is increased in two ways. ae 
first is through latent heat release, primarily Dy convec- 
momo cin the warm sector of tne storm. Second, large 
IMitial air-sea temperature differences result in a strong 
surface heat flux which acts to maintain the low-level 
horizontal temperature gradient. In addition, the sensible 
heat added within the PBL is rapidly mixed to mid-levels 
mimougnearyeconvective adjustment. Rapid intensification 
Ceo rOleenvemoaroclini€s conversion of tnis AE to KE. 
The amount of AE directly generated by latent heat release 
in the warm sector, where eddy efficiency Factors are 
Mesitives, 1S about an order of magnitude less than that 
Peoaiiced py eddy heat transport. Overall, the inclusion of 
diabatic effects increased the kinetic energy of the polar 
low by over an order of magnitude. 

Mii vedatumCnnIsS Generated In the dilabatic simulation 
of tne polar low by Kelvin-Helmholtz instability. Maximum 
Vac smOumcitoministabtlity are Tound in tne trough asso- 


Ciated with the polar low. While tnis instability does not 


Ze 


appear to be significant to the overall development of the 
polar low, it may be significant in tne formation of meso- 
scale features near the front as suggested Dy Bennetts and 
Hoskins (1979) and observed by Locatelli et al. (1982). 

The simulated polar low, with diabatic effeets Tncioece 
ies f ounemres resemble a simulated synoptic-scale extra- 
tropical cyclone (wavenumber 8), as observed by Locatelli et 
al. (1982). Similarities exist Init ec rnicmom ven tempera- 
ture, precipitation and surface pressure. The effects of 
diabatic processes on wavenumber 8 is found to be similar to 
that.for the polar Mowe 

The results of this sttidy confirm the major yee = 
that the Pacific polar Jow is a Daroclimirc adi sam oaneee 
However , <Gilidibiatike “Cage tales), specifically latent heat release 
Tne the Warm Sect ore cainicimet) spi e sc Mer mol rcla initial @iressed 
temperature difference of at ledstmos Ky are Neeessaryenon 
the polar low to achieve observed deepening rates. This 1s 
consistent with the speculation of Sardie and Warner (1983) 
that moist bdroclinity 1s responsible jor the qrowc neo) an 
Pacific polar low. 

Research should be continued to further our under- 
standing and prediction capabilitles og. ene pal cmmnor. 
Specifically, our understanding of convective and planetary 
boundary layer physics must continue to improve since they 
are so crucial to the development oj Syne monic ailleoteeeene ee 


can be accomplished through Observ de oma les Gade ee ere 


Lie) 


studies are needed to determine the precise role of Kelvin- 
Pommnelez Instability mMmeaqr the cold front and its possible 
ele nin the formation Of rainbands. Also, this study does 
not address either the Atlantic polar low or the continental 
popecee OW. Studies sSimliar@to this one, wsing= the approe- 
Memece Tmltial Conditions, dre necessary to obtain a better 


understanding of these systems. 
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APPENDIX A: THE FORECAS ia Dian 


A. DYNAMICS 
The forecast model uses the primitive equations in the a 


coordinate system after Phillips (1957), 


Bae E 
oh ie (A-1) 
where 0°< o ( 1. The Vel@tatieismdnicr 
a(mu) _  a(muu)  _(muv) _ 3(muc) an 
on 3x ay 30 POTTY Sear 
A-2) 
_ ae fogee | 
ae tr - TK, ue (Vou) 
d(tmv) . _ d(tuv) a(tmvv) — d(mvo) oT 
at 3x - oy 30 pies KY ay 
34 ae (A-3) 
=a TF + TK, ve CO) 
o(mT) _ 3 (muT ) 9(mvT) 3(1oT) RT1o Cale 
ae 3X : ay - Ye) > “cit aaa Wm) 
A-4) 
RT Ze TT ( 
ct cC [ VetWdo + TeV eed = ieee eel 1: C 
p*o p 
a(mq) . _ a(muq) _ a(mvq) _ a(moq) agar 
at aX dy ele} m 
(A-5) 





= = oi V7 (mV )do (A-6) 


0 
Nee Se we 
SuGe = J Dee ide - Vao-(rV) (A-8) 
0 


iemvept Camo Ewe enGe Of tne Grid 1S SnhOWnM In Fig. A=1. 


All variables are defined at the middle of each layer, 


except nO, Witeme:S Gdrrled de stheunterface between the 


layers. 


The sigma levels are defined in Table A-1l. 


Taormina «6©Model levels. 
Level Number O 
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24 
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The horizontal structure of the Gelduis Shown “In F1d: 


A-2. This is defined as scheme C by Arakawa and Lamb 


leg 7 


interval 


Tor comnicmU—componenc 15) deétTined one-half grid 


in the x-direction from the mass points and 


ical 


J=0 
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Vertical structure of model grid with distribu- 


tion of forecast variables. 


eZ 


pole, =2 3 HOnrnzZOneccimrrstri1bDUcCTOn Of Mass and momentum 
HOnnesmam tne forecast grid. 
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the v-component is defined one-half grid interval in the 
y-direction from the mass points. [he nid open came 
in x and has free-slip walls on the y boundaries. 

The model uses the split-explicit time integration 
scheme developed by Madala (1981). First, the equations are 
integrated over a time step governed by the meteorological 
and slowest gravity modes. The solution is then adjusted 
for the linear terms which govern the fastest gravity modes. 
The momentum and thermodynamic equations are solved with 
centered time and space differencing, while the moisture 
equation differencing is forward in time and upstream in 
space. Fourth order advection is used for the momentum 
and temperature. Second order advection is used for the 
mol stuire-. 

Fourth order diffusion is used in all the ohare keh keen ne 
equations except for nm. Second Order diffi us von) se camer 
the first interior row. This diffusion 1s done On gner sre 
levels, which can lead to spurious sources and sinks of 
energy when applied to the mass and moisture fields. fo 
overcome this problem, the di fiUsTOnme iowa po neds bomen 
tions from the standard atmosphere for the temperature and 


moisture. The difftlsion Coenii cient si smecmou team. 


Bia eas on lst ihberi omeeroy 
Ky = <1.6 AX On 2nAd-4t hime enero (A-9) 


1.25AX On 21) Other. Wiaie mi ore ccm 
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mee hxe 1s tne grid Ssodcing In meters. High frequency time 
Semmiations dre controled bY dpolying the Robert (1966) 


mien ht) cen with tne smoothing coefficient set to 0.15. 


eer at SICS 
Peecianetdry  sbounddry Layer 
The surface fluxes are calculated in the manner 
fecertoed bY Deardortf (1972). First, the stability of the 
PBL is calculated using the bulk Richadson number, 
Rig = ghin vs) (A-10) 


B 


Um 


lew ala ene vertical potential temperature, suoscript 


Ss denotes surface values and subscript m denotes mean PBL 


momes. For stable or neutral conditions (Rip ae i, Ene 
friction coefficient is 
Rig 
Cy = OF (1 - RT)» (A-11) 
iC 
ame the meat transfer coefficient 15s 
Rip 
Cy = Con (l-5-—) (A-12} 
a 1G 
where the neutral values are 
0.025h =a 
C = [ —~ In (— 3) + 8.4 ] (A-13) 
uN iM Z 


Z R Shel Size a - 

Ce, Ly ioe ) + 7.3 ] (A-14) 
and 

Ric = ene (A-15) 
Rip was constrained to be no larger than 90% of Rig. The 
surface roughness, Lo is a specified constant over land and 

2 
LZ. = Uso ue (A-16) 


over water. The value of Lo Over water is constrained to be 
no less. than.2.435 © POs mm (5 Ugg co qimmn wii oem ONa a 
NEPRF). 


For unstable comd? tions hj CO eee si mieetesons 


B 
GQeCTT1elente) > 


C, = e = 25 @€xp | U5 2p ote sy (A-17) 
and the he@t transfer coetiie venus: 
- - | - | -l,-] 
C . = (C, + on 1° oe ) : (A-18) 
where 
SEG py ak t) = 258 (A-19) 
and C., a midi sate Gaye ieee (A-13) and (janes 


Gesmeceively. Once C and C, apemenown, themsurface stress 


momecdtculated from 


— 


ME eee (A-20) 


gee ere Virtual temperature flux 15 
Ue (ee eS aa (A-21) 


Deardorff (1972) proposed extrapolation equations to 
determine the mean PBL Gucci ee Sama hiyis cdl) lead ¢€O Sem nouUS 
Meoeems 11 tne PBL does not extend at least through the 
Boevomemodel layer. 10 avoid extrapolations when the PBL 
is shallow, the PBL is constrained to always be at least as 
deep ACE hic lowest model layer (about 40 mb). As long as 
the lowest layer is reasonably shallow, any errors intro- 
duced should be relatively small since such a shallow PBL 
will be presumably stable and the surface fluxes should be 
relatively small. 

The PBL is assumed to always be well-mixed. Each 
time the model physics routines are called (this can be any 
Mmamegem mulclpleswo, the dynamic time step), u, v, @ and q at 
all model levels within the PBL are replaced by average 


Zoesethrougnout the PBL. That is, 


OpBL 
ag 


2 S22) 


Carre 


Vo 7 


where @ Can De Us y, Semone. By making these assumptions, 
the values found in tne lowest layer of the model are always 
the medn quantities. 1 nwtiem, cee 

The height of the PBL in unstapte cone ee donc ei 
modeled after Stull (1976). The prognostic equation for 
the PBL depth (aP) is 


5 | (A-23) 
Pea meme UD )) W + W 


arc 

(a) Poy ie , 
Term (a) defines the effect of large-scale motion on the 
PBL. Term (b) is the entrainment velocity given by 


O.1(w' . 
eee ot 


e AG 


vm 
Oy 
m Z 3 
+ gd, A8, [0.05u, |W | a C00) Cave) i (A-24) 
m 
where 
2 
Uy 
Ace = (A=72 oe 


Mele 16 the inversion strength at thew: seam dq, 
is the depth of the layer of negative heat flux at the top 
of the P Biee 


WAG an 
d) = St ___ , (A-26) 
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The entrainment velocity is found by iterating (A-24) and 
P20 Ihe TirSt term On the right-hand side of 

(A-24) describes the buoyant production of turbulent 
Sti cechiergy, while the second term describes the mechan- 


Pomme rodvelvon Of energy. Stull's equation for W also 


contained a term for energy loss due to internal gravity 
waves. This term is neglected here. Term C in (A-23) 
accounts for cloud-induced subsidence. The computation of 
Mieceeerm 1S described later in this Appendix. 

The prognostic equation for the PBL depth in stable 
Gemd1ciOns tollows the development of Nrewstadt and 


Tennekes (1981). 





(A-27) 


The equilibrium height is given by 


a 2 
0.156, fsina cosa |V)| 


Ne = g[d8,7dt] os) 
over land when [de /dt | SO eee enr= Sand 


| bu 


ns F flv [sina 


(A-29) 


over land or water when |dg /dt| < 0.1°K hr-~! after Brost 


and Wyngaard (1978). The cross-isobar angle is arbitrarily 


pe 


specified as 10° at all times. The tiiwemiegciiopeeseee 7 


is Given Dy 


0.75(6, -96_) 
5 = oes (A-30) 
dey Gt 
fo) 
when (8, - coe > 1°K, and by Khakimoues ( 197Geeeeai atime 
s = O87 (A-31) 
when (6) 90) < digak 


2. Ree 1 Omnitele 


A modified version of the Kuo (1965) cumulus 
parameterization is used. This modified version links the 
convection to the PBL. First) the molstune conv eng enc yim 
is calculated as 


My = Ve(q mW) (1-dpe, ) 7 Pen he (A=3 2 


\, 
g 
Here, Mz 1S the moisture convergence occurring solely in the 
PBL. Convection can occur when M+ > 4.0 X 1075 gm m-és-l 
and when the equivalent potential temperature decreases with 
height from the PBL to the first model layer above the PBL. 
The final constraint 1s that Convect 1 on we animseoc Cun ene 
LCL 1s above the PBL. In addition, 1f Convection 1seoeeune 
ring, the cloud-induced subsidence terme £ermG eae 


(A-23) is computed as 
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Neo = Mh cy > Opp, )/2At (A-33) 


Micmenueis the time step tor the physics routine. This 
feeowto Ggrive tne PBL Naltway down to tne LCL in the absence 
Smmoener effects. 


Cloud temperature (T ara mpxTnd fa) Os (qo) are 


mecumed £0 be those=of the moist adiabat which passes 
mg@emon the LCL. the Dartitioning of the neating and 
moistening follows Kuo (1965). A value b is defined so 
that a fraction DM, is used to moisten the environment and 


(1-b)M, condenses and falls instantaneously as rain. The 


Mawes Or b 1s defined as 
O 

D = Og (1-IRH)do . (A-34) 
9 


where IRH is the integrated relative humidity (0 < IRH < 1} 


eee, ene Geoth of the cloud, 





al 
[RH = — | (-L)de (A-35) 
CmamnnGee 


Nomen cet mu GeG]OTLatiyOn 15 assumed to occur 


lemons deieau1on CONd]t10NS O¢CCUr at a Grid. point. 


14] 


Moisture is condensed until saturation 1s reached. All 
condensed moisture falls Into the next layer wniecnetseenem 
checked for supersaturation. Precipitation occurs only when 
the air is saturated t7om the €l clld “comriewe ound. 
3. ery Convect | vieeeciiisiemieme 

Ory convective adjustment is applied whenever the 
lapse rate EASE the dry adiabatic lapse rate. Tempera- 
tures in the unstable layer are modified such that the total 
potential energy is conserved. All unstable layers are 


adjusted siibgmeect to tne Condit ome 


meres (A-36) 


eX) 
J0 





and 


(A=398 


fey 


[ 
(Sie dom cae 
ot 


where AT 1s the temperature change required to reduce the 
lapse rate to neutral. These equations are solved as a 


system of linear equations. 


aly Incorporation of Physica poeeoacs 


Physical processes tend to evolve relatively slowly 


and are not restricted as severely by time step limitations 
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as are the dynamics. Therefore, the physical parameteriza- 
PPCnMCwOt ley eT Ocd tO Oeteaetilated at some multiple of the 
iames. time SteD. sin this study, the non-convective 
heating and the cumulus parameterization are applied with 

a 2400s time step. The heating and moistening rates are 
distributed uniformly over all time steps until these 
moumcines are called ag@din. fhe PBL processes and the dry 


convective adjustment are called every dynamic step. 
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APPENDIX B: ENERGY BUDGE) SOA iii 


The energy equations presented here were developed Dy 
Robertson and Smith (1983). Several boundary terms do not 
appear in this study because of the closed domain. The 


equations use the Tollowingdefiirewouc. 
ee | (Mice (eae 


is the zonal average, 
(--) = 4 a (B-2) 
is the deviation from the zonal mean or the eddy component, 


and 


Ca ela )dy (B-3) 


<|— 
of 
< 


is the average in the y-direction. The mass-integrated 


kinetic energy 15 


| V-V 
K = leak (—z) dp . (B-4) 


and the mass-integrated available potential energy is 


T 
A = Cm | [ET]dp (B-5) 
O 
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mere wei S "une Efficiency tector, defined as, 
————— (B-6) 
ive reference pressure, oe fret nerareal mean pressure on an 


isentropic surface. The zonal and eddy components of the 


Kinetic energy are 


TT ee 
V 
K2 = 7. id (3-7) 
O ; 
and 
TT 
7] y'.V' 
KE = _ [—5-—]dp (B-8) 
0 


respectively. The zonal and eddy components of the 


available potential energy are, 


Ce: 
.2 mT [ET]dp (B-9) 
@) 
and 
Cc AY 
cen Lif [leet dl|cayo Rt 0 } 
0 


respectively. 
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The budget equations are: 


where, 


Gye 
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moO, Sima 0 








Lw'a'jdp 


CA “GA Geeare 


Gee peers 
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{ Lw aldp 
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TT 
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(B-16) 


(B-l 
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DZ (8-20) 
0 
in 
C ae 
seme a CE'OQ'Jdp (8.2 2ie 
‘9 
TT 
C aaa 
gv = a [EQidp (B22) 
QO 
a oT 
PAE = fhe All ees Oech) 
oy eta 
PAL = 3 be or! Bene) 
View 
=) Lp 23 Sy : 
PKE : [ 5 ie (8-25) 
7 ] VW. OT 
PKR i ae (B-26) 
Peni Olmeowl iGersea description Of these terms: CA is the 


conversion of energy between AZ and AE; CE is the conversion 
eeeeondyvecectween) Ae and KE by Vertical redistribution of 

masS, which is often referred to as baroclinic development; 
CK is the conversion of energy between KZ and KE, eee aS 
MeseierenenrG~ea CGedas tne barotropic term; CZ 1S the conver- 


Sion of energy between KZ and AZ by vertical redistribution 
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of mass; DE (DZ) 3S the dissipation OF Kemer ce 
friction; GE (GZ) is the generation of AmegZ eee tc ciemee 
processes; PAE (PAZ) is the change in AE (AZ) due to changes 
in the surface pressure and PKE (PKZ) is the change in KE 


due to cnanges in the surface pressure. 
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APPENDIX C: NONLINEAR SIMULATIONS OF WAVENUMBER 8 


iiemolrgese On Lilsmaunpendix 18 to show the results of 
iembanear, dadilabatlie anid diddatic Simulations of wavenumber 
Perec cegemes is) jetatoebe compared to the results for 
ci lan Simulations ofewavenumber 20 ‘presented in Chapter 
IV. The minimum central pressures for the adiabatic (G8A) 
and diabatic (G8D) experiments are given in Fig. C-1l. The 
central pressure in G8A falls from 1003 mb at 72 h to 968 mb 
ieee O Mewich rising pressure afterwards. in G8D, the 
central pressure falls from 999 mb at 72 h to 944 mb at 
oe ne Ties deepening rate is similar for both systems from 
72-108 h. However, after this time, G8D exhibits a rapid 
decrease in central pressure while the central pressure in 
Gepe StdresetOrrise slowly. tne reasons for this pressure 
tendency will now be exatnined. 

icweec—terve Oressure Tields for 60, 84, 108 and 132 h 
ZieeconweGdmachnare shown im Figs. C-2 and C-3, respective- 
ly. At 60 h, both forecasts exhibit an open trough pattern 
ee Liew Dresciure C-jmmgb lower in GaO. In addition, the 
mEPougi —extenasmrurthewm tO themsoutheast in G8D witn the 
Pemmacten iO: sd seCcOnd Center sinmeinis trougn. The extension 
Meelis or OUugGwmPsa due EO Cemvectilon in the warm sector of 
Siewaiaclnuodniicenmm oy of 1 GOA NaS attained a closed circu- 


lation with a central pressure of 996 mb, whereas G8D has 
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deepened to 991 mb. By 108 h, both systems have intensified 
further with G8A near the mature stage. At 132 h, GSA 
remains at nearly the same central pressure while G8D has 
continued to deepen and has reached the mature stage with a 
Seema pressure of 947 mp. 

Cross-sections of height perturbations through G8A and 
fowemoiercs fh are given ig Fig. C-4. Both systems exhibit 
maximum perturbations at the surface and in the upper levels 
with a westward tilt with height. The upper level height 
feeeerae ion 1S Stronger an Geb than in G8A, which shows 
miduedidbatic efrects tend to strengthen the upper level 
Secearion.  iINis is alee foUnd to be the case for the 
cmecpemtows (Cit 1g.04-5 em IKxe =G2UA, which is strongly 
damped in the vertical, the upper level perturbation in G8A 
is more than half the magnitude of the Surface perturbation. 
Tnis again illustrates that shorter waves tend to be shallow 
Gineiemances as TOUNd itm Chapter [I1. 

Cross-sections of vertical motion, temperature and the 
PagiewmEaxiS TOreGon dnidmogumat S4 9h are shown in Fig. C-5. 
Tne adiabatic case exhibits a warm perturbation extending 
from the surface to 400 mb. The center of the warmest air 
Ti oaneNdiiewltn NeigGhteprom tne Surface to 500 mb, and 
westward ae 500 mb. A center of upward vertical motion 
OCCURS JUSt Gast of tnhemeGrface low in the Warmest air, 
which implies an increase of KE over a deep layer due to 


maim eenbecer Clm@elidt ion. Tats 15 17 Contrast to 
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the polar low, G20A (Fig. 4-2), which exhibits very little 
tilt to the warm air perturbation from enews iin de oe ome 
with a strong westward tilt above 7/50 mb. 

Significant differences exist between the cross-sections 
of GBA and G8D. First, theveenter of upward Vertycaleieeao 
and the warmest air in G8D are found about 7/50 km eastuor 
their respective positions in GSA. This 1S a "1 ost) gee 
convection taking place ahead of the surface low as noted 
above. Second, in G8D, a narrow band of upward motion 
within a weak warm air perturbation is found in the cold air 
about 1000 km behind the position of the Surface lOws Ine 
Surface reflection of this feature is the weak pressure 
trough in the northweSt=,low 1ngowmet S29 h (Rig. (C— im 
This pressure trough is a weak polar low and its formation 
is related to an interaction between convection and the PBL. 
Forecast PBL depths in the cold air behind the trough are 
relatively large because the maximum air-sea temperature 
differences are in this area. Some convection (about 5 mimn/ 
12 h) is forecast behind the suriace slow in the cold aimee 
This results in a decrease in the PBL height at the convec- 
tive points due to cloud=indwcedmaubesd cnicem ba. A= 33) 
the points where the PBL is deep, relatively strong PBL 
winds exist due to the mixing of the winds throughout the 
PBL. At the convective points, where the PBL is shallow, 
relatively light winds exist siiieemeciic: Stimondemewindseed Were 


are not being mixed into the PBU tines omen cee 7 ce Or 


ya 


convergence in the PBL on the western edge of convective 
teas waren fTurtner enhances Convective forcing of this 
yin sii We Tals eopceesse 1s ecapable of producing smail- 
Secens cold dirtmdiseureances in this model, it is felt that 
Milemis mare Of a numerically tmduced feature as opposed to 
Mm crea ly realistic situation. ~Certainly, not al] 
convection drives the PBL halfway down to the LCL every half 
hour as is assumed here. However, large variations in PBL 
depth over a small area in the presence of strong vertical 
wind shear snould be considered as a possible mechanism to 
trigger polar low development. 

ile meowd| IecOlmvceLIve and Staplesprecipitation accumu- 
lated in G8D during the 12 h period ending at 103 h are 
ween FiGesc-o. ADOUL 2ZOumm of precipitation is forecast 
to the north of the disturbance with a maximum of 38 mm to 
PiemeoOurMNedSt "Ofecne TOweealomg themwcold front. Only about 
40% of the northern maximum and 2/3 of the southern maximum 
In precipitation is convective. Therefore, the percentage 
of convective precipitation is smaller in G8D than in the 
polar low, G20D. 

The effect of the diabatic processes can be seen in 
Plossl -7, sana ©-6 in which the 850, 700 and 500 mb height 
and temperature fields are given for G8A and G8D, respec- 


tively. As with the polar low in Chapter IV, a noticeable 
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warming occurs ahead of the developing wave due to the 
release of latent heat in this area, whicn reduces tne 
Cimect Of adiabatic coolimg through ascent. 

The surface sensible and latent heat fluxes associated 
WeeeeGol dt WO@en are shown 1n Fig. C-9. Maximum values of 
the fluxes are found behind tne developing wave in the 
ompcie dir, dsemoted with the polar low (cf, Figg 4-11). 
However, unlike the polar low, downward (negative) fluxes 
Pweranodd OT Goleman ENe warm Sector. This results from a 
large low-level temperature perturbation (Fig. C-8a) in G8D 
which reverses the sign of the air-sea temperature differ- 
ence in the warm sector. 

KE budgets for G8A and G8D, integrated over the forecast 
domain, are presented in 6 h intervals wn Tables C-1 and 
Co-2erespectcively. ihe maximum KE attained by G8A, 6.7 x 
ei Scmabemt “O¢ hess Ehametne Maximum attained by 
G80, 11-1 x 10° Jm-2. This is in sharp contrast to the 
Pooes in Cmameer [TV that the diabatic effects ied to an 
order of magnitude increase in Wel or G20D. The reason why 
G8D does not exhibit a similar increase in magnitude over 
G8A can be seen from Fig. C-10, which gives the time-height 
eee eros Or GsheanduGol. For the polar low 
diseteccdmlieemaoter [Ys the adiabatic run exhibits a maxi- 
mum of KE near the surface with negligible amounts above 
600 mb. However, G8A exhibits a maximum KE of 5x 104 J m-2 


near 350 mb with a secondary maximum of 3 x 104 J m-2 near 
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900 mb. Therefore, the reason G20D reali zecmsU cae apc 
increase of KE when dij@ibatic efm@eiets anemia lieder eect 
G20A is a very shallow system and the domain-integrated KE 
is quite small. WttQhidvabatic effeets. the pola meow 
extends through a ldmge part of the tropospmere and yreld. 
a much larger domain-integrated KE. G8A already extends 
through the troposphere and the only otner way to increase 
the domain-integrated KE is to increase the circulation of 
the system. This o@eurs in GSD and results 1h kK Sema cen 
9x 104 J m-2 at 35Qgmb and 6 x PG moog O0mnioe 

The baroclinic development term, CE, 1s the largest term 
in the budgets for DoOuhGGA Sana Gav tChimoughmes2 fh. Oe ig 
the developmental stages prior (tomeaIeeh te is l=l/25touz 
times larger in G8D and imethe 1 argesgaeterm in the Ke 
budgets. Time-height cross-sections of CE are given in 
Fig. C-ll. Both foreedasts exhibi taiaximum=econVersion ae 
to KE at /00 mb, whiten is near the level op mda. imu ere cc 
motion. Large differences exist between the vertical struc- 
ture of CE for G8A and G20A. G20A exhibits maximum CE near 
900 mb with negative values above /00 mb. This is a result 
of the warm air perturbation tilting sharply westward above 
7/50 mb. However, themwarm dir perce ia bon 1 ac Gee 
eastward with height up to 500 mb during the development 


stage. This results i posh Give valUucsmo; CHlup toes0eem. 
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in G8A. The vertical structure Of Ce Mise eee ee 
between the diabatic simulations of G20 (cf, Fig. 4-13) 
and G8. 

The horizontal Darotrvepien term, HWChKe 1 Suecemenneonn commen 
both forecasts. Maximum Values of HCK are found ay time | 
mature stages of G8A and GOD] 8] itme-netqnc cross-sec uvoncmom 
HCK are given in Fig. C-12. The largest values of HOCK are 
found in the upper levels associated with the jet. Also in 
G8D, positive values of HCK are found below 600 mb from 96- 
114 nh. The reason for this “3 the southeastward extension 
of the developing trougn forced Dy COnVeCclioOn as Seen need 
Surface pressure field for GSD (Fig. €-3). nis pat tein 
implies a positive correlation between u' and v’ and hence 
a4 conversion. of Ke Ceuta Ae yaic 

The vertical component of CK, VCKe 1s posit i Veuenine wan 
out both forecastsy@and therefore represents asourcemonmiass 
As with the polar low Forecasts, VCK 1S amaxi mum G= 12a 
after CE reacnes its maximum value near tne mature stage of 
the cyclone. This should be expected since u'’ and w' are 
largest then. The fact that VCK attains its largest values 
at or just after maximum storm Intensity suggests senate an 
process is a symptom of the baroclinic disturbance rather 
than a cause; as 78 also the case foOyeHGK. [herenenc woe 
G20 and G8 are considered tO bed istunmmancesmanich forms 


through conventional Daroci imiG sims clo ie ee ie — iene in 
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cross-sections of VCK for G3A and GSO aye Oi cae. 
C-13. Both disturbances exhibit maximum values ™oT Vilieeo, 
approximately tne same magnitude between 850-900 mb and near 
400 mb. 

Frictional dissi@ieati0n, DE, aes (as curs iO pen commmnne 
both forecasts. DBE 15 the largest termine ween oo cicero 
Once the disturbance reaches itS.maximum intensity. Tnis is 
in contrast to G20D, which exhibited larger values of HCK 
than DE during the decaying stage of the polar low. How- 
ever, the trends here are similar Since for GeoAy 9UE = 1sece 
times larger than HCK at 132 heWwoai bei or Gc Ue DE cee 
times as large as HCK at the same time. Therefore, 
barotropic damping is larger relative to DE when diabatic 
effects are included. 

KE budgets integrated over the forecast domain are 
presented in 6 h intervals for G8A and G8D in Tables C-3 and 
C-4, respectively. As was the case for the polar low, G8D 
exhibits more AE than G8A. The maximum values for G8A, 2.2 
x 105 J m-2, and fae Gop. 4.3 x Wt eo cena tle some 
Time-height cross-sections of AE for G8A and G8D are given 
in Fig. C-14. Altnotigh Gepmexh 1 bts lah Glen ey artes mon er ancre 
all levels, the most significant increase over G8A occurs in 
the mid-levels from 700-500 mb. AE values for G8D in this 
layer are three times larger from 120-132 h when the central 
pressure fell 7 mb in G8D and rose 2 mb in G8A. It was 


noted earlier that the largest differences in CE between G8A 
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Mimocimeccur Inmrne 700-500 mp layer from 120-132 h (Fig. 
C-ll). Tnis is at the same time and location in which the 
large increases in AE are observed, which suggests that the 
disturbance uses the additional AE for baroclinic develop- 
ment of the storm. 

iemocdoman anmenseoy  diabDatic processes, GE, attains 
fiaeaimineyaldes trom !i42138 Hh. This 1S consistent with the 
large increases in AE noted during this time. The time- 
WenummemecrOSS-section of GE for G8D is given in Fig. C-15. 
The largest values from 114-138 h are found in the 700- 

500 mb layer. This is attributed to the latent heat release 
primarily in the warm sector ahead of the low as discussed 
above. Negative values of GE in the lowest 100 mb are the 
result of positive (negative) sensible heat fluxes occurring 
iimewe =oolde warmoair behind {ahead of } the lowe: Weis 
Decor imet GemGesemp |\@smtnat tor G20) (Fig. 4-17). 

The horizontal eddy heat meamspont, HCA, 1S PGSHmuve for 
Dem can waicmookomnaien indicates thateit is a source of AE. 
During the developing stages of G8A and G8D, HCA is the 
largest term in the AE budget. Values of HCA are up to two 
times larger in G8D, which suggests that the diabatic 
processes have altered tne overall temperature distriodution 
Cmarniemilsuiigwadanieesso @S CO Increase the horizontal eddy 
heat transport. Jhis was shown in Chapter [V to play a 
Significant role in the intensification of the polar low. 


Time-height cross-sections of HCA for G8A and G8D are given 
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Tae ride G=!6. Although GeA stiggests a maximum near 450 mb 
and another near 300 mb, both forecasts exhibit large values 
throughout the 900-300 mb layer. A similar pattern is found 
for G2Z0D except over a slightiy shallower layer. 

Noticeable differences exist in the vertical eddy heat 
Eta OT ¢ oe VCA, between G8A and G8D. In the adiabatic 
run, VCA is positive during the developing stage up to 102 
amen indiCatesmemat VCA 15 a source @f AE. After 102 4h, 
VCA 1S negative and tends to damp tne wave. However, renee 
Of VCA are generally 4-10 times less than the largest term, 
ewe Gob. VChe1s Mositive through 90 h. While G8D 1s 
deepening rapidly from 96-144 h, VCA is negative (except at 
P20 @) wand 1s Very smaliaein magnitude. Even though VCA is 
(ec cdi ines tinismmmne. the values are not as large in 
magnitude as in GSA, and therefore more AE is available in 
G8D. Time-height cross-sections of VCA for G8A and G8D are 
given in Fig. C-17. Generally, negative values are found 
above 500 mb. Below 500 mb, G8A is dominated by positive 
values prior to 108 h and negative values after 108 h. The 
Same tendency exists in G8D except that positive values 
persist at all times below 800 mb. 

ieee er TOTmOi AL by VCATin GSA and GSD during the 
growth stage is different from that found in the polar low. 
In GAGpy Generavgion Gi Ee Bc nrougneyCAsduUring»the developing 
Stage was found to be approximately the same magnitude as 


term HCA. The reason for tnis appears to be the area in 
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which the low develops. VCA is comprised of the four 





terms + a wal E/Cy ; b ) wEta/C, ; Cc). (Gh dE /tm amassed 





(ou'€) see Term ais generally an order of magnitude 
less than tne otner terms, so VCA is a balance between term 
DS which is typically less than zero, and terms c and d, 
which are typically positive. The relative strength of eacn 
term governs the sign of VCA. In G20D, term b tends to De 
relatively small Simee tims dustin bance; onic] Sumo cu mman 
farther to the north) than in Go -and therefore @ isecmaerers 
This reduces the magnitude of term b and VCA is positive. 

The results of the G8 simulations presented here illus- 
trate tne following poOunts2 a oe dfjapatic processes serve 
to intensify an ext#atropicd!, SVymopeee-ceale diStintcmecumn 
a similar manner to the Pacific polapslow.. Soeecitieaen = 
the diabatic processes serve to alter the temperature struc- 
ture of bDotn G8 and G20 such that Gddyeaeat transpor eae 
therefore values of AE, are increased. Tne direct genera- 
tion of AE by diabatic processes (Gb |) Mieeound to beac 
least 5 times smaller than the eddy heat transport in botn 
G8 andjG20. Second?) patterns of SUrtideemmnressiire, pweenen— 
tation, wind and temperature for a Pacific polar low are 
found to resemble those of synoptic-scale oceanic 
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